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1.  PRODUCTION  AND  DETECTION  OF  H4* 

1.1  Ha*  REVIEW 

The  four-center  H4  system  has  been  the  subject  of  theoretical  studies  since  the 
late  sixties,  primarily  due  to  interest  in  the  relationship  of  the  ground  state 
hypersurface  to  the  H2  +  D2  ->  2HD  exchange  reaction. 1 * 2 , 3  Investigations  into 
the  nature  of  the  excited  state(s)  of  such  system  have  been  limited  in  number 
even  though  SCF-MO-CI  calculations  in  1976  by  Goddard  and  Csizmadia  indicated 
the  existence  of  a  bound  "excimer  state". ^  Recently,  however,  Nicolaides  et 
al.,  using  an  extension  of  their  modeling  work  on  rare  gas  hydrides^  have 
predicted  the  existence  of  a  bound  singlet  excited  state  of  H4  and  have 
extended  this  analogy  to  the  ^2)3  and  ^2)5  clusters. 

The  theoretical  studies  of  Refs.  6  and  7  predict  that  the  bound  state  of  H4 
arises  from  reaction  of  molecular  hydrogen  in  its  ground  state  with  H2  that  has 
been  excited  to  the  B(1EU+)  state. 

HjtfB1  V)  +  ^(X1  £g+)  +  M  -»  H4*  +  M  .  (1) 

Further,  the  stable  excited  state  is  predicted  to  be  produced  only  when  the 
B^TU+  state  precursor  has  sufficient  "ionic"  character  for  reaction  with  H2(X). 
This  is  predicted  to  occur  for  v  >  3  in  the  upper  state. 

Quenching  studies  of  HD(B,v)  performed  by  Atkins,  Moore,  and  coworkers®-^  show 
an  unusually  large  quenching  cross  section  (79A^)  for  H2.  This  had  led  to 
further  speculation  that  formation  of  a  complex  is  involved.  The  H2(B)  state 
in  v'  =  3  has  sufficient  energy  to  dissociate  itself  as  well  as  an  H2  collision 
partner: 

H2(B)  +  H2  ->  4H  AE  =  -2.70  eV  . 
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Hence,  quenching  of  H2(B)  by  H2  may  be  highly  reactive  with  H-atoms  as  the 
products. 

The  calculations  of  Nicolaides  and  coworkers^-?  show  that  H4*  is  bound  by 

3.1  eV  with  respect  to  dissociation  to  form  H2(B)  and  H2.  However,  they  also 
have  shown  that  the  complex  is  bound  by  1.86  eV  with  respect  to  dissociation 
along  the  H-  +  H3+  coordinate.  Hence,  photolysis  of  H4*  may  produce  H3+  +  H~, 
H3  +  H,  or  H2(B)  +  H2-  A  careful  selection  of  photolysis  wavelengths  may 
result  in  selective  excitation  into  any  of  these  channels.  In  our  experimental 
effort  we  will  be  using  a  range  of  photolysis  wavelengths  in  order  to  excite 
the  channel  leading  to  H2(B)  +  H2. 

Little  has  been  published  concerning  the  structure  and  bond  frequencies  of  the 
H4*  excimer.  The  structure  is  predicted  to  be  trigonal  pyramidal,  comprising  a 
H3+  complex  in  an  equilateral  trigonal  geometry  of  sides  0.899A  and  and  H~ 
center  above  the  H3+  triangle  at  a  bond  distance  of  2.01A.  The  electron 
density  of  the  structure  is  higher  in  the  center  than  along  the  bonds.6  Though 
bond  vibrational  frequencies  for  H4*  have  not  been  calculated,  the  structure  of 
the  H3+  center  is  quite  similar  to  that  calculated  for  free  H3+.  The 
calculations  of  Carney  and  Porter^  have  shown  symmetric  stretch  and  bending 
vibrational  frequencies  for  the  ion  to  be  3185  and  2516  cm-*.  Hence,  in  both 
laser  photodissociation  and  passive  emission  studies  of  H4*  progressions  in 
these  frequencies  may  be  observed.  Passive  emission  from  the  H4*(^A')  origin 
may  appear  as  progressions  in  the  H2<X,v)  system  at  frequencies  shifted  by  the 

3.1  eV  binding  energy  of  H4*  with  respect  to  the  H2(B-X)  Lyman  system.  Thus  a 
careful  study  of  both  the  absorption  and  emission  spectrum  of  the  excited 
mixture  is  required  to  detect  H4*  formed  from  the  reaction  of  82(B)  with  H2. 

1.2  DESIGN  OF  EXPERIMENTAL  APPARATUS 


The  goal  of  the  design  phase  of  the  experimental  program  was  to  produce  a 
device  which  would  enable  us  to  produce  H2(B,v')  in  a  supersonic  jet,  provide  a 
sufficient  number  of  collisions  with  H2  to  produce  H4*,  and  probe  for  its 
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existence  and  molecular  properties.  We  chose  to  produce  H2(B,v)  in  the  jet  by 
energy  transfer  from  Ar*  resonance  states  (1048A,  1066A  lines)  in  Ar/H2 
mixtures  via  discharge  production  of  Ar*  in  the  nozzle  expansion  region. 

Detection  of  H4*  may  be  via  observation  of  passive  emission  from  the  radiating 
excimer  or  via  photodissociation  of  H4*  to  produce  H2<B),  which  may  be  observed 
in  emission. 

1.2.1  Nozzle  Source 

The  primary  reason  for  using  a  supersonic  molecular  beam  to  produce  H4*  is  to 
provide  a  locally  high  density  but  optically  thin  source  of  H2  for  direct 
excitation  of  the  H2(B)  state  as  well  as  collisional  production  and 
stabilization  of  H4*.  Three  primary  criteria  must  be  met  by  the  source: 

1.  There  must  be  sufficient  collisions  in  the  expansion  beyond  the  nozzle 
throat  for  collisional  quenching  of  H2<B)  to  occur  at  a  rate  at  least 
comparable  to  that  for  spontaneous  radiation. 

2.  When  energy  transfer  from  Ar*  to  H2  is  employed  there  must  be  sufficient 
collisions  beyond  the  nozzle  threat  for  Ar*  to  be  quenched  by  H2. 

3.  Formation  of  stabilized  H4*  is  a  three-body  process.  Hence,  sufficient 
collisions  must  occur  in  the  post-excitation  region  between  the  colliding 
H2<B)  -  H2  complex  and  a  third  body  (H2  or  Ar). 

A  further  practical  constraint  is  that  the  total  gas  flux  from  the  system  must 
be  handled  by  a  2000  ls“l  diffusion  pump.  This  was  primarily  a  cost  constraint. 

These  constraints  are  met  using  pulsed  valves  of  the  types  reviewed  by  Gentry.12 
These  valves  typically  have  gas  pulse  durations  of  100  to  500  us  at  repetition 
rates  up  to  35  Hz.  These  duty  factors  are  ideally  suited  to  experiments 
employing  pulsed  lasers.  Due  to  the  reduced  gas  load  realized  by  using  these 
sources,  smaller  pumping  systems  may  be  employed  to  achieve  desired  background 
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pressures.  Typically,  the  desired  background  pressure  is  set  such  that  the 
molecular  mean  free  path  is  comparable  to  dimensions  of  the  background  chamber. 
In  our  case  this  is  approximately  lO-^  Torr. 

Gas  flow  in  freely  expanding  supersonic  jets  has  been  treated  extensively  by 
Anderson  and  and  Fenn.-*-^  These  sources  are  characterized  by  a  rapid 

drop  in  number  density,  a  decrease  in  kinetic  temperature,  and  an  increase  in 
beam  translational  energy  downstream  of  the  nozzle  throat.  These  properties 
are  typically  expressed  in  terms  of  the  Mach  number  of  the  flow,  M,  and  the 
distance  downstream  of  the  nozzle  throat  in  nozzle  diameters,  (L/D).  The 
number  density,  N,  the  temperature,  T,  and  the  velocity,  V,  are  given  by: 


N 


V  = 


[l  +  (Y  -  1)  M2/2 j 

[l  +  (r  -  1)  M2 / 2 ] 
1/2 

M 


l/(l-r) 


-1 


(2) 

(3) 

f4) 


where  y  is  the  heat  capacity  ratio  of  the  gas  (Cp/Cv)  and  the  subscript  zero 
denotes  the  nozzle  plenum  or  "stagnation"  conditions.  The  value  of  M  must  be 
determined  numerically  and  has  been  tabulated  by  Anderson^  as  a  function  of 
(L/D)  for  various  values  of  y.  The  rapid  drop  in  number  density  and 
temperature  in  these  expansions  is  illustrated  in  Fig.  1. 

We  will  use  expressions  (2)  to  (4)  to  define  the  nozzle  operating  conditions 
subject  to  the  collisional  quenching  constraints  imposed  by  our  optical  pumping 
methods.  There  is  some  question  concerning  the  degree  to  which  rotational 
degrees  of  freedom  couple  to  translation  in  H2>  We  have  assumed  they  couple 
well  in  our  calculations  (y  =  7/5)  as  opposed  to  y  =  5/3.  In  the  case  of  Ar* 
quenching  by  H2  to  form  H2<B)  a  quenching  rate  coefficient  of  (1.9  ±  0.5)  x 
IQ-10  cm3  molecule-*  s~!  has  been  measured  by  Boxall,  Chapman,  and  Wayne. 
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Figure  1.  Supersonic  nozzle  reduced  number  density  and  temperature  profiles 
for  monatomic  (y  =  5/3)  and  diatomic  (y  ■  7/5)  gases. 

The  radiative  rates  for  the  Ar*  resonance  lines  are  given  by  Viese  et  al.^7  as 
5.1  x  108s_l  (1048A)  and  1.19  x  108s~l  (1066A).  The  condition  for  efficient 
quenching  is 


kQ[H2]  >  k 


Ar 

rad 


(5) 


If  the  faster  radiative  rate  is  employed,  then  ve  require  [H2I  >  2.7  x 
1018  molecules  cm-8  tor  efficient  quenching  or  Ar*.  For  H2(B)  quenching  by  H2 
the  results  of  Fink,  Atkins,  and  Moore8  give  a  super-gas-kinetic  quenching  rate 
coefficient  of  1.1  x  lO-^  cm8  molecule-^  s~l  at  300  K.  The  radiative  rate  for 
the  B-state  has  been  measured  by  Hessler^-8  to  be  1.25  x  lO^s"^-.  Hence  H2  con¬ 
centrations  in  excess  of  1.1  x  lO^-8  molecules  cm-8  are  required  to  quench 
H2(B).  Thus,  at  H2  densities  sufficient  to  quench  Ar*  radiation  H2(B)  will  also 


5 


be  quenched.  If  the  radiative  rate  for  the  Ar*  1066A  line  is  used  in  expres¬ 
sion  (5)  the  H2  quenching  density  obtained  is  6.3  10^  molecule  cm-^.  A  useful 
practical  H;  ’uenching  density  is  thus  approximately  1  x  10*8  molecule  cm-^. 


Average  mass  flow  through  the  nozzle  is  a  function  of  the  nozzle  stagnation 
pressure,  cross  sectional  area,  flow  velocity  at  the  throat,  and  valve  duty 
cycle.  The  pumping  speed  required  to  achieve  a  given  background  pressure  is 
simply  written  as 

Kp  ■  S  £  Tp  °2  vo  "o'nb  <6> 


where  f  is  the  valve  pulsing  frequency  (10s~^),  Tp  is  the  pulse  duration  (10_/|s), 
D  is  the  nozzle  throat  diameter  (cm),  Vq  is  the  beam  velocity  at  the  throat 
(M  =  1,  V  =  1.2  x  105  cm  s-*),  Nq  is  the  stagnation  number  density,  and  Ng  is 
the  desired  background  number  density  (10-^  Torr  =  3.3  x  10^  molecule  cm~3). 

The  values  of  these  parameters,  where  provided,  are  typical  for  the  nozzle 
source  under  consideration.  Figure  2  shows  the  pumping  speed  required  to 
achieve  10"^  Torr  operating  pressure  for  four  stagnation  pressures  (1,  2,  5,  and 
10  atm)  as  a  function  of  nozzle  throat  diameter.  The  horizontal  line  on  the 
plot  defines  the  pumping  speed  of  the  Varian  VHS-6  diffusion  pump  system  appro¬ 
priate  to  our  system.  The  intersection  of  the  diffusion  pump  line  with  the 
pumping  speed  carves  define  the  appropriate  operating  nozzle  diameters  for  each 
pressure.  Using  these  values  of  Nq  and  D  we  can  use  the  curve  in  Fig.  1  to 
determine  the  distance  from  the  nozzle  throat  at  which  the  [ H2 ]  is  sufficient  to 
quench  both  ^(B)  and  Ar*.  This  is  plotted  in  Fig.  3  and  shows  that  the  dis¬ 
tance  is  approximately  0.9  mm  and  is  independent  of  the  combinations  of  Nq  and  D 
which  require  the  maximum  pumping  speed.  Thus  0.9  mm  is  the  maximum  distance 
from  the  nozzle  throat  at  which  H2(B)  can  be  excited  and  efficiently  quenched. 

Ve  have  chosen  to  employ  a  5  atm  nozzle  stagnation  pressure  with  a  250  um 
diameter  throat  to  continue  the  design  calculations  since,  in  our  previous 
experience,  this  combination  has  been  shown  to  be  optimum  for  mechanical 
performance  of  the  nozzle  under  consideration.  The  H2  expansion  centerline 


6 


PUMPING  SPEED  (6s-i) 


STAGNATION  PRESSURE  (atm) 

A-4860 

Figure  3.  H2<B)  quenching  distance  as  a  function  of  nozzle  stagnation  pressure. 

number  densities,  kinetic  temperature,  and  beam  velocities  are  shown  for  these 
conditions  as  a  function  of  distance  from  the  nozzle  throat  in  Figs.  4  to  6. 

The  H2(B)  and  H4*  collision  rates  may  be  calculated  and  integrated  from  these 
plots.  The  molecular  collision  rate  is  simply  given  from  kinetic  theory  as: 


Z  =  l2  iucN  (7) 

where  a  is  the  collision  cross  section  and  c  is  the  average  molecular  velocity. 


Hence,  the  collision  frequency  is  a  function  of  number  density  and  pressure, 
both  of  which  are  decreasing  with  distance  from  the  nozzle  throat.  The  total 
number  of  collisions  a  molecule  undergoes  is  simply  the  time  integral  of  the 
collision  frequency,  which  is  decreasing  with  distance  from  the  throat.  In 
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Figure  4.  H2<B)  nozzle  expansion,  H2  centerline  number  density 


Figure  5.  H2<B)  nozzle  expansion  kinetic  temperature. 


9 


EXPANSION  DISTANCE  (mm)  ^ 

Figure  6.  H2<B)  nozzle  expansion  beam  velocity. 

divided  by  the  beam  velocity,  which  is  also  changing  with  distance.  This  inte¬ 
gration  and  the  other  nozzle  operating  parameters  are  shown  in  Table  1.  The 
integration  shows  that  the  H4*  complex  may  undergo  300  collisions  if  formed 
0.9  mm  from  the  nozzle  throat.  This  assumes  a  relaxation  cross  section  identi¬ 
cal  to  the  H2(B)  relaxation  by  H2 ( 7 9A^ ) .  Stabilization  of  the  H4*is  essen¬ 
tially  a  two  step  process.  This  first  step  is  the  collisional  stabilization  of 
the  orbiting  H2<B)  -  H2  complex  during  its  lifetime.  The  second  step  involves 
quenching  of  the  highly  excited  H4*  to  lower  vibrational  levels  of  the  complex. 
The  first  step  is  the  most  critical  because  it  requires  that  a  collision  occur 
in  the  10"^-®  to  10“12s  during  which  the  excited  complex  exists.  This  requires 
very  high  collision  rates.  Our  estimate  is  that  H4  formation  efficiencies  in 
the  range  of  10~2  to  10_<4  per  H2<B)  molecule  excited  may  be  expected.  This 
would  result  in  approximately  5  x  10^  to  5  x  10^  H4  molecules  being  produced 
for  each  laser  pulse.  Detection  sensitivities  are  difficult  to  estimate  due  to 
lack  of  any  information  about  the  H4*  photodissociation  cross  section.  If  this 
cross  section  is  in  the  range  from  10”^®  to  10-^®  cm^  then  only  modest 
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TABLE  1.  Nozzle  performance  characteristics. 
Total  collisions  =  306;  L  >  0.9  mm. 
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TOTAL  COLLISIONS  *  306 
L  >  0.9  mm 


11 


photodissociation  fluxes  are  required  to  produce  ten  counts  per  laser  pulse  at 
the  signal  averager.  Hence  all  three  of  the  design  criteria  fc-c  the  nozzle 
have  been  satisfied  for  this  configuration.  The  H2(B)  collision  frequency  is 
shown  in  Fig.  7  as  a  function  of  expansion  distance.  The  H2(B)  radiative  rate 
is  given  by  the  dashed  line.  This  confirms  that  excitation  at  distances  les^ 
than  0.9  mm  are  required  for  H2<B)  quenching  to  dominate  over  radiative  decay. 

The  molecular  beam  source  we  have  chosen  for  our  application  is  manufactured  by 
Newport  Research  Corporation  (Model  BV-100/BV-100D) .  This  is  the  commercial 
version  of  a  design  by  Adams  and  coworkers^  and  is  shown  schematically  in 
Fig.  8.  The  valve  operates  at  stagnation  pressures  up  to  10  atm  with  a  pulse 
duration  as  low  as  100  us  and  a  repetition  rate  of  50  Hz. 


EXPANSION  DISTANCE  (mm) 

AU608 

Figure  7.  H2 ( B )  collision  frequency  as  a  function  of  distance  from  the  nozzle 
throat . 
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Figure  8.  Cutaway  view  of  solenoid-actuated  pulsed  beam  source  (Adams  et  al., 
1981).  (1)  Actuator  assembly  with  dual  solenoids  and  ferrite  cores. 

(2)  Nozzle  plate.  (3)  Mounting  ring. 

1.2.2  Discharge  Excitation 

Discharge  excitation  of  Ar*  in  Ar/H2  mixtures  was  accomplished  using  the 
throat-discharge  modifications  to  the  NRC  pulsed  beam  valve  as  described  by 
Grant  and  coworkers^O  and  further  modified  by  Physical  Sciences  Inc.  (PSI). 

These  modifications  are  shown  in  Fig.  9.  The  basic  concept  involves  the  passage 
of  the  expanding  gas  through  a  region  between  two  biased  plates.  The  two 
plates  are  separated  by  a  ceramic  insulator.  Grant  and  coworkers  experienced  a 
problem  with  feedback  from  the  discharge  affecting  the  pulsing  electronics  for 
the  nozzle.  To  eliminate  this  problem,  a  second  ceramic  insulator  was  placed 
between  the  discharge  plate  and  the  nozzle.  This  isolates  the  discharge  ground 
from  the  nozzle  and  eliminates  most  of  the  feedback  problems.  At  higher 
discharge  voltages,  feedback  in  the  form  of  double  pulsing  of  the  nozzle  was 
still  a  problem.  Discharge  voltage  was  always  increased  while  monitoring  the 
pulsing  signal  to  the  nozzle.  If  a  second  pulse  was  observed,  then  the  voltage 
was  decreased  until  the  nozzle  operated  normally. 
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When  the  gas  nozzle  is  closed,  the  pressure  is  sufficiently  low  that  no 
discharge  occurs  between  the  plates.  As  the  valve  opens  and  gas  fills  the  void 
between  the  plates,  initially  spark  breakdown  occurs  between  the  plates  and 
then,  as  the  pressure  increases,  a  transition  to  glow  discharge  conditions 
occurs.  Each  of  these  conditions  has  been  treated  extensively  by  Cobine^  for 
the  case  of  breakdown  between  plane  parallel  electrodes.  Our  discharge  con¬ 
figuration  is  different,  but  it  is  useful  to  carry  out  these  calculations  to 
roughly  parameterize  these  conditions. 


Sparking  voltages  for  H2  as  a  function  of  pressure  times  distance  are  shown  in 
Fig.  10.  For  our  gap  distance  of  3  mm,  it  can  be  seen  that  only  modest  voltages 
are  required  to  cause  the  gas  to  arc  early  in  the  expansion.  Extrapolation  of 
the  H2  curve  in  Fig.  10  to  background  conditions  (pd  *  3  x  10"^  Torr-cm)  shows 
that  very  large  voltages  are  required  to  cause  arcing  under  these  conditions. 


VOLTS 


pd  (torr-m) 

A-4613 

Figure  10.  Paschen  curves  for  H2,  CO2,  and  air. 

The  transition  to  glow  discharge  operation  occurs  when  the  gas  pressure  exceeds 
the  voltage  required  for  arcing.  At  this  point,  the  gap  voltage  must  be 
sufficient  to  sustain  a  glow  discharge.  The  condition  for  sustaining  a  glow 
discharge  is  given  by: 

y(e°d-  1)  =  1  (9) 

where  a.  is  the  coefficient  for  field-intensified  ionization  by  electrons,  y  is 
the  probability  for  electron  production  at  the  cathode  from  positive  ion 
bombardment,  and  d  is  the  gap  spacing  in  cm.  The  coefficient  a  is  a  function 
of  pressure  and  gap  voltage.  A  plot  of  oc/P  (cm-*  Torr"l)  as  a  function  of  E/P 
(V  cm-*  Torr~l)  is  shown  in  Fig.  11  for  H2,  the  dominant  component  in  the  Ar/H2 
expansion.  The  value  of  y  is  relatively  insensitive  to  voltage  for  slow  moving 
ions  and  has  a  value  of  0.053  for  Ni  electrodes . 21  Expression  (9)  can  be 
solved  for  this  value  of  y  and  d  =  0.3  cm  to  give  a  value  of  a  =  9.96  eta"!  to 
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A-4814 


Figure  11.  Coefficient  of  field-intensified  emission  for  H2. 

sustain  the  glow  discharge.  Using  the  curve  of  Fig.  11  and  the  range  of 
stagnation  pressures  possible  for  use  in  the  valve,  the  plot  of  minimum  bias 
voltage  versus  gap  pressure  shown  in  Fig.  12  is  obtained.  In  the  discharge 
excitation  experiments,  the  discharge  is  conducted  in  the  throat  region  and 
stagnation  pressures  as  low  as  0.4  atm  may  be  employed  to  achieve  the  required 
collision  frequency  at  the  point  of  H2(B)  formation.  Under  these  conditions,  a 
bias  voltage  of  1.5  kV  is  required  to  sustain  the  discharge. 

The  glow  discharge  current  is  limited  to  a  few  hundred  microamps  by  placing  a 
3  MS  ballast  resistor  in  series  with  the  spark  gap  and  the  power  supply.  The 
two  parallel  discharge  plates  have  a  calculated  capacitance  of  4.7  pf.  A 
0.0047  uf  capacitor  was  added  in  parallel  to  increase  the  total  capacitance  of 
the  discharge  circuit.  The  ability  to  monitor  the  discharge  was  accomplished 
by  adding  a  1  k2  resistor  in  series  with  the  voltage.  We  could  then  determine 
when  the  discharge  fired  reproducibly  under  various  conditions. 
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P“<atm> 

Figure  12.  Minimum  gap  voltage  required  for  sustained  glow  discharge  in  H2  as 
a  function  of  nozzle  stagnation  pressure. 


1.2.3  Discharge  Characterization 

The  gas  valve  is  fired  through  an  internally  generated  pulse.  The  discharge 
then  self-fires  as  the  gas  pressure  increases  between  the  two  plates.  Under 
these  conditions,  the  discharge  pulse  has  a  2  us  duration.  This  is  consider¬ 
ably  shorter  than  the  several  hundred  microseconds  the  valve  is  open.  There  is 
also  50  us  of  firing  jitter  between  the  pulse  that  fires  the  gas  valve  and  the 
firing  of  the  discharge.  Attempts  to  reduce  the  jitter  by  adjusting  discharge 
voltage,  gas  pressure  and  repetition  rate  were  unsuccessful.  To  eliminate  the 
jitter  from  our  detection  electronics,  the  signal  from  the  discharge  monitor 
was  used  as  a  trigger  source. 

Since  the  discharge  pulse  is  shorter  than  the  valve  gas  pulse,  the  question 
arises  as  to  when  in  the  gas  pulse  the  discharge  is  firing.  The  gas  valve  is  a 
solenoid  valve  requiring  50  us  to  open  fully  and  reach  full  gas  pressure  at  the 
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throat.  If  the  discharge  fires  before  the  valve  fully  opens,  then  the 
concentration  of  H2  could  be  considerably  less  than  what  is  required  for 
quenching  of  H2<B). 

A  discharge  voltage  of  -380V  was  needed  to  achieve  breakdown.  The  discharge 
was  operated  at  700  to  800V  for  a  pure  H2  discharge.  A  slightly  lower  voltage 
was  needed  for  f^/Ar  mixtures  to  prevent  feedback  to  the  pulsed  valve.  The 
discharge  proved  to  be  fairly  insensitive  to  pressure  which  could  indicate  the 
discharge  fires  before  the  valve  fully  opens.  A  pressure  of  400  Torr  was  used 
to  assure  a  H2  concentration  large  enough  to  quench  H2(B). 

The  question  arises  as  to  whether  we  ever  achieved  glow  discharge  or  only 
created  a  spark  discharge.  We  calculated  a  voltage  of  1.5  kV  was  required  for 
glow  discharge  at  400  Torr  of  hydrogen.  However  we  could  not  operate  at  this 
voltage  because  of  feedback  problems.  If  the  valve  was  not  opening  completely 
before  discharge,  the  actual  pressure  would  be  less  than  400  Torr  and  a  glow 
discharge  formed.  The  pressure  insensitivity  of  the  discharge  would  indicate 
that  this  might  be  true  for  our  operating  conditions. 

I. 2.4  Detection  of  H4* 

Detection  of  H4*  was  attempted  via  observation  of  passive  emission  from  the 
excimer  and  via  photodissociation  to  form  H2(B).  Passive  emission  could  be 
observed  using  either  a  vacuum  monochromator  system  (Acton  Research  Model 
VM-502)  or  an  ^-purged  broadband  radiometer  (130  to  190  nm).  Nicolaides^-? 
and  coworkers  have  calculated  the  binding  energy  of  H4*  to  be  3.1  eV  with 
respect  to  H2(B)  and  H2(X)  products.  Since  the  energy  of  H2(B)  is  about 

II. 7  eV,  the  excimer  might  be  expected  to  emit  at  energies  less  than  8.6  eV  or 
145  nm.  These  wavelengths  may  be  easily  observed  with  an  ^-purged  mono¬ 
chromator  system.  A  Deslandres  table  for  the  H2(B-X)  transition  is  shown  in 
Table  2.  Correlation  of  the  wavelengths  in  this  table  with  the  Franck-Condon 
factors  for  the  transition  calculated  by  Spindler^  shows  that  emission  from 
H2(B)  excited  at  the  nozzle  throat  will  also  be  observable  using  this  detection 
system.  We  can  use  this  emission  as  a  diagnostic  to  determine  when  we  are  in 
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TABLE  2.  Deslandres  table  giving  wavelengths  in  nm  for  H2(B-X)  transitions. 


0 

1 

2 

3 

H2(B,m 

4 

3 

4 

7 

8 

0 

110.340 

109.244 

107.733 

106.327 

104.973 

103.489 

102.470 

101.314 

100.217 

1 

114.222 

114.448 

112.313 

111.249 

109.748 

108.343 

107.034 

105.773 

104.579 

2 

121.778 

119.354 

118.041 

114.329 

114.711 

113.179 

111.729 

110.356 

109.034 

3 

127.517 

125.408 

123.425 

121.533 

119.789 

118.120 

114.541 

113.047 

113.435 

4 

133.417 

131.111 

128.944 

124.904 

124.981 

123.145 

121.449 

119.328 

118.297 

5 

137.448 

134.931 

134.349 

132.349 

130.258 

128.287 

124.427 

124.671 

123.015 

4 

145.347 

142.824 

140.259 

137.349 

133.582 

133.448 

131.434 

129.340 

127.752 

7 

151.717 

148.742 

143.940 

143.332 

140.902 

138.398 

134.430 

134.387 

132.465 

3 

157.323 

134.404 

131.402 

148.791 

144. 153 

143.474 

141.347 

139.154 

137.096 

H21X,V) 

9 

143.791 

140.329 

157.101 

134.084 

151.257 

148.406 

146.115 

143.775 

141.577 

10 

149.308 

145.803 

142.333 

139.133 

154. 

153.296 

150.448 

148.162 

145.323 

11 

174.339 

170.900 

147.237 

143.322 

lid. 631 

157. a44 

154.844 

152.219 

149.757 

12 

179.433 

173.478 

171.418 

’>8.024 

144.449 

161.531 

138.593 

153.340 

153.260 

13 

183.724 

179.380 

173.34. 

171.598 

148.100 

144.831 

161.773 

138.909 

154.228 

14 

184.941 

182.445 

178.277 

174.401 

170.709 

167.416 

164.262 

161.310 

158.548 

15 

189.117 

184.317 

180.233 

174.294 

172.403 

169.159 

163.939 

162.928 

160.111 

14 

190.103 

185.434 

181.130 

177.130 

175.  .23 

149.948 

144.698 

163.660 

160.817 

17 

189.783 

185.151 

180.859 

176.872 

173.158 

169.692 

144.452 

143.422 

160.588 

18 

188.085 

183.334 

179.314 

173.396 

171.743 

168.332 

145.144 

162.161 

159.370 

resonance  with  ^(B-X)  transitions  as  well.  A  lens  system  was  used  to  image 
the  H2  beam  centerline  on  the  monochromator  slits  for  improved  signal  strength. 

A  boxcar  signal  averager  was  employed  to  record  the  signal  from  the 
monochromator  system. 

1.2.5  Overall  Chamber  Design 

The  vacuum  chamber  was  comprised  of  an  8  in.  diameter  four-way  cross  (Fig.  13) 
mounted  on  a  diffusion  pump  stack.  The  diameter  of  the  cross  is  matched  to  the 
diffusion  pump  throat  to  minimize  conductance  losses  through  the  chamber.  The 
pumping  system  consists  of  an  8  in.  diffusion  pump  (Varian  VHS-6),  a  liquid  N2 
cooled  cryotrap  (Varian  316),  and  an  aluminum  slide  valve.  The  total  through¬ 
put  of  the  system  is  2040  Is-*.  All  joints  in  the  system  are  either  ASA  o-ring 
or  knife-edge  flanges.  The  molecular  beam  valve  is  directed  horizontally  in 
the  chamber.  Eight  ports  located  symmetrically  around  the  nozzle  can  be  used 
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Figure  13.  Schematic  diagram  of  H4*  apparatus. 

for  laser  access  and  detection  of  H2(B)  for  H4  fluorescence  vith  a  solar  blind 
photomultiplier.  To  resolve  the  fluorescence,  a  vacuum  ultraviolet  monochroma¬ 
tor  is  attached  to  the  end  of  the  vacuum  chamber  opposite  the  pulsed  nozzle.  A 
CaF2  lens  is  used  to  image  the  throat  of  the  nozzle  onto  the  slit  of  the  mono¬ 
chromator.  For  our  experiments,  it  was  sufficient  to  purge  the  monochromator 
vith  nitrogen. 

1 . 3  RESULTS 


1.3.1  H?(B)  Spectrum 


A  spectrum  of  the  emission  from  the  discharge  region  vas  obtained  using  the 
vacuum  ultraviolet  monochromator  and  recorded  on  a  computer.  A  typical 
spectrum  for  a  pure  H2  discharge  is  given  in  Fig.  14.  The  spectrum  vas  taken 
vith  10A  resolution.  A  spectral  fit  to  the  experimental  spectrum  is  then 
generated  by  computer  in  three  steps.  First,  a  very  high  resolution  spectrum 
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Figure  14.  ^(B)  emission  spectrum.  Thin  line  represents  measured  spectrum 

and  thick  line  is  the  spectral  fit. 


covering  the  wavelength  range  and  vibrational  levels  needed  is  generated  by 
calculating  the  emission  intensities  and  line  positions  for  a  given  molecular 
system.  Second,  vibrational  basis  sets  are  generated  by  passing  an  appropriate 
slit  function  over  the  high  resolution  spectrum.  Finally,  after  the  data  is 
corrected  for  instrument  response  and  an  average  background  intensity  is  sub¬ 
tracted,  a  least  squares  fit  of  the  experimental  data  to  the  basis  functions  is 
calculated.  The  best  fit  to  the  data  is  given  in  Fig.  14.  A  rotational  temp¬ 
erature  of  1000  K  vas  used  to  obtain  the  fit.  This  is  reasonable  considering 
the  hot  plasma  nature  of  the  discharge.  Since  the  resolution  of  the  spectrum 
is  so  low,  individual  vibrational  transitions  are  hard  to  distinguish.  A 
vibrational  distribution  can  be  determined  from  the  spectral  fit  and  is  given 
in  Fig.  15.  The  relative  electron  excitation  efficiency  for  producing  excited 
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Figure  15.  H2<B)  vibrational  distribution  as  determined  from  calculated 

spectral  fit  to  experimental  spectrum. 

H2(B)  is  governed  by  the  Franck-Condon  factors  which  connect  H2(X)  with  ^(B) 
in  the  absence  of  additional  feed  mechanisms.  The  resulting  vibrational  dis¬ 
tribution  produced  by  Franck-Condon  type  excitation  is  also  given  in  Fig.  16. 
The  two  differ,  with  the  discharge  peaking  at  much  lower  vibrations.  This  non- 
Franck-Condon  distribution  could  be  due  to  vibrational  relaxation  in  the  valve 
throat  or  an  additional  production  channel  of  low  vibrations.  It  cannot  be 
accounted  for  by  excitation  of  H2<X,v)  at  a  vibrational  temperature  less  than 
or  equal  to  the  electron  temperature. 

1.3.2  H?/Ar  Mixtures 

Additional  H2<B)  should  be  produced  by  co-discharging  H2  +  Ar. 


Ar  +  e" 

->  Ar*  +  e" 

(10) 

H2  +  Ar* 

•*  H2(B)  +  Ar 

(ID 
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Figure  16.  Variation  of  H2(B)  intensity  with  variation  in  Ar/H2  mixture. 

Our  calculations  have  shown  that  Ar*  will  be  quenched  at  similar  H2  con¬ 
centrations  needed  to  quench  ^(B).  Monitoring  the  H2(B)  fluorescence  signal 
while  adding  Ar  to  the  discharge  mixture,  we  observed  a  decrease  in  the  H2(B) 
signal.  The  variation  of  H2<B)  intensity  with  varying  amounts  of  Ar  is  given 
in  Fig.  16.  The  signal  is  largest  for  the  pure  H2  discharge,  then  decreases  as 
Ar  is  added  and  finally  plateaus  at  large  Ar  concentration.  No  evidence  for 
Ar*  energy  transfer  to  H2  to  form  additional  H2(B)  was  observed.  This  effect 
can  be  explained  by  considering  that  we  are  co-discharging  H2  +  Ar.  In  this 
system,  there  is  a  competition  for  the  available  electrons.  Any  electrons  used 
to  form  Ar*  are  no  longer  available  to  produce  H2(B).  This  electron  depletion 
by  argon  is  much  more  efficient  then  energy  transfer  from  Ar*  to  H2. 

1.3.3  Atomic  Hydrogen 

Spectra  of  the  discharge  were  also  recorded  in  the  visible.  The  spectrum  given 
in  Fig.  17  shows  several  peaks.  This  can  easily  be  assigned  to  atomic  hydrogen 
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Figure  17.  H-Atom  emission  spectrum  from  discharge  source. 


transitions.  Since  an  absolute  calibration  of  our  detection  system  has  not 
been  done,  no  comparison  can  be  made  of  the  ratio  of  atomic  hydrogen  to  H2(B) 
formed.  The  quenching  rate  coefficient  of  H2(B)  by  atomic  hydrogen  is  unknown 
but  may  be  an  additional  depletion  channel  of  H2(B). 

1.3.4  Ha*  Production 

Formation  of  stabilized  H4*  is  a  three-body  process.  Sufficient  collisions 
must  occur  in  the  post-excitation  region  between  the  colliding  H2(B)-H2  complex 
and  a  third  body.  By  forming  82(B)  at  the  throat  of  the  gas  valve,  we  should 
have  sufficient  collisions.  However,  if  the  discharge  is  firing  before  the 
valve  opens  fully,  sufficient  collisions  may  not  occur.  After  confirming  that 
82(8)  was  being  formed  in  our  discharge  source,  ;he  next  step  was  to  try  to 
detect  the  formation  of  H4*. 
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Detection  of  H4*  may  be  via  observation  of  passive  emission  from  the  radiating 
excimer  or  via  photodissociation  of  H4*  to  produce  H2<B),  which  may  be  observed 
in  emission.  Emission  was  detected  with  a  broadband  radiometer  (130  to  190  nm) 
or  a  monochromator  (190  to  540  nm).  All  emission  observed  from  the  discharge 
could  be  assigned  to  H2(B)  or  atomic  hydrogen;  no  passive  emission  due  to  H4* 
was  observed. 

A  nitrogen  laser  was  used  for  photolysis  of  H4*  at  337  nm.  The  photolysis 
laser  was  delayed  with  respect  to  the  discharge  source.  Delays  of  0  to  30  us 
were  used,  but  no  H2(B)  emission  was  observed. 

There  is  a  second  channel  possible  for  H4*  photolysis: 

H4*  -»  H3+  +  H' 

Though  we  cannot  currently  detect  ions,  this  channel  does  need  to  be  examined. 
However,  a  different  source  of  H2(B)  would  be  needed  since  our  discharge  source 
produces  ions.  A  cleaner  source,  such  as  laser  excitation  of  H2(B),  would  be 
ideal . 

1.3.5  Detection  Limits 

Since  we  were  unable  to  observe  any  signal  from  H4*,  it  is  important  to  estab¬ 
lish  detection  limits  on  our  system.  First  we  need  to  calculate  the  [ H2 ( B ) 1 
formed  in  our  source.  This  is  determined  by  the  electron  density  of  our  dis¬ 
charge  source,  the  excitation  rate  coefficient,  and  the  concentration  of  H2. 
Using  the  current  measured  by  our  current  monitor,  the  electron  flux  for  the 
discharge  area  can  be  calculated  to  be  1.25  x  10^-®  e-  cm-^  s-*.  Ajello  et 
al.23  have  measured  the  emission  cross  section  of  H2(B)  by  electron  impact  over 
an  energy  range  from  threshold  to  400  eV.  The  threshold  for  this  system  is 
20  eV.  In  our  discharge  the  average  energy  is  assumed  to  be  2  eV.  The  frac¬ 
tion  of  electrons  above  threshold  is  then  calculated  assuming  a  BoltEmann  dis¬ 
tribution.  The  excitation  rate  is  then  the  product  of  the  cross  section  and 
the  fraction  of  the  electrons  above  threshold.  For  our  operating  conditions, 
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we  calculate  an  excitation  rate  coefficient  of  0.106  e“/s.  The  concentration 
of  H2  is  1.29  x  1019  molecules  cm-^.  The  production  of  ^(B)  is  then 
1.4  x  1018  molecules  cm"8  s-l.  The  discharge  occurs  for  2  us  so 
[H2(B)1  =  3  x  10l2  molecules  cm“3. 

The  next  problem  is  to  determine  how  much  H4*  will  be  formed  from  H2(B).  The 
relevant  equations  for  this  process  are: 


H2(X)  +  e-  H2(B, v) 

Jl  -  »I 

(12) 

H2(B)  +  h2  -»  h4* 

<j)k2 

(13) 

H4*  +  M  -»  2H2  +  M 

k-3 

(14) 

H4*  +  M  -»  H4  +  M 

k3 

(15) 

H4*  ->  H2(B)  +  h2 

k4 

(16) 

H2(B)  ->  H2  +  h\> 

a5 

(17) 

The  H4*  lifetime  is  very  short  so  a  steady-state  solution  is  appropriate: 


d[H4] 

3T 


=  J1[H2] 


'  ♦  k2[H2]  ' 

'  k3[H2] 

tal82l  *  *SJ 

Zk3lH2]  +  k4 

(18) 


The  assumption  has  been  made  that  k_3  =  k3  where  k3  is  the  lifetime  of  H4*. 

One  can  then  estimate  the  H4*  formation  efficiency  per  H2(B)  and,  using 
t H2 ( B ) ]  =  3  x  10l2  molecules  cm“3,  calculate  an  estimated  [ H4* J  for  a  predicted 
H4*  lifetime.  The  results  are  given  in  Table  3.  The  dissociation  cross 
section  of  H4*  is  uncertain  but  a  reasonable  estimate  would  be  10“ 19  cm2#  For 
our  laser  at  337  nm  we  could  dissociate  0.02  percent  of  the  H4*.  The  actual 
photodissociated  volume  is  0.2  cm^.  If  one  assumes  that  all  dissociated  H4 
forms  H2<B)  then,  for  an  H4*  lifetime  of  10"10s,  dissociation  would  form 
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TABLE  3.  Production  of  H4*. 


k3  (s'1) 

H4*  Quenching 
Efficiency 

[H4*]  Produced 
(molecules  cm-3) 

Estimated 

h2(B)^h4* 

Branching  Ratio 

1010 

10-2 

3  x  1010 

0.02 

1011 

10-3 

3  x  109 

0.20 

1012 

10-4 

3  x  108 

>  1 

2.5  x  10?  H2(B)  molecules  in  0.2  cm3.  Since  we  are  now  in  the  collisionless 
regime  of  the  gas  jet,  only  the  radiation  channel  is  available.  The  radiative 
lifetime  of  H2(B)  is  1.25  x  109  s~3  and  the  laser  pulse  duration  is  10  ns. 

This  gives  a  total  of  3  x  10**  photons.  The  number  of  photons  can  be  converted 
to  a  current  by  considering  the  collection  efficiency,  transmission  of  the  lens 
and  filter,  gain  of  the  photomultiplier  and  quantum  efficiency  of  the  photo¬ 
cathode.  For  our  system  3  x  10**  photons  would  give  a  signal  current  of 
6.4  x  10~2A.  The  signal-to-noise  limit  of  the  system  is  1  x  10~**A.  Therefore, 
if  the  lifetime  of  H4*  is  shorter  than  5  x  10~*-2s,  we  would  not  be  able  to 
detect  it.  In  other  terms,  if  the  excited  state  lifetime  of  H4*  is  100  ps,  the 
efficiency  for  producing  H4*  from  quenching  of  H2(B)  is  less  than  0.02. 

In  this  argument  we  have  made  the  major  assumption  that  all  dissociating  H4* 
forms  H2<B).  If  this  were  not  true,  our  detection  limit  of  H4*  would  be  even 
higher,  as  would  our  yield  estimate. 

1.3.6  Summary 

The  discharge  nozzle  source  does  produce  H2<B)  at  the  throat  of  the  nozzle. 
However,  the  H2<B)  is  formed  in  several  vibrational  states  and  rotationally  hot 
(-  1000  K).  No  passive  emission  from  H4*  was  observed.  An  attempt  was  made  to 
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photolyze  any  H/4*  present-  to  form  Ht(B).  No  additional  Hq(B)  fluorescence  due 
to  the  photolysis  laser  was  observed. 

1 . 4  QUENCHING  OF  H?(B) 

The  goal  of  this  program  was  to  selectively  excite  H2(B)  in  a  supersonic  jet  in 
such  a  way  as  to  promote  the  formation  of  the  H4*  excimer  from  collision  of 
H2(B)  with  H2.  It  is  important  that  the  H2(B)  be  quenched  by  H2  before  it  can 
radiate.  The  radiative  rate  for  the  B-state  has  been  measured  by  Hessler^  to 
be  1.25  x  10^  s~*.  Some  limited  data  exists  on  the  quenching  rate  coefficient 
of  H2 ( B )  by  HD.  Fink,  Atkins,  and  Moore®  give  a  super-gas-kinetic  quenching 
rate  coefficient  of  1.1  x  10^  cm®  molecule"*  s~*  at  300  K. 

In  this  study  we  have  measured  the  quenching  rate  coefficient  for  H2(B)  v'=3 
and  v'=4.  For  each  vibrational  level  we  have  measured  the  quenching  rate 
coefficient  for  several  different  rotational  levels. 

In  the  calculations  of  Nicolaides  et  al.,®  they  predict  that  the  stable  excited 
state  of  H4  is  produced  only  when  the  H2  B-state  precursor  has  sufficient 
"ionic"  character  for  reaction  with  H2(X).  This  is  predicted  to  occur  for  v  >  3 
in  the  upper  state.  With  laser  excitation  of  H2,  we  can  produce  H2(B)  in  a 
single  vibrational  level.  However,  there  is  the  possibility  of  vibrational 
relaxation.  Therefore  we  have  examined  whether  any  vibrational  relaxation 
occurs  during  quenching. 

1.4.1  Experimental  Apparatus 

H2(B)  is  formed  by  three-photon  laser  excitation  of  H2.  The  fluorescence  from 
H2<B)  is  then  measured  as  a  function  of  hydrogen  pressure.  The  apparatus  used 
for  the  quenching  measurements  is  shown  in  Fig.  18.  Laser  light  of  the  proper 
wavelength  is  generated  by  a  Nd:YAG  laser  pumping  a  tunable  dye  laser,  which  is 
then  frequency  doubled.  The  laser  light  is  then  focused  into  a  vacuum  cell  con¬ 
taining  a  known  pressure  of  hydrogen.  To  resolve  the  H2(B)  fluorescence,  a  vac¬ 
uum  monochromator  with  a  solar-blind  photomultiplier  (EMR  Photoelectric  Model 
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Figure  18.  Schematic  of  H2<B)  quenching  experimental  apparatus. 

5420  is  used  to  detect  the  signal.  The  photomultiplier  is  attached  directly 
to  the  cell  when  total  fluorescence  detection  is  desired.  The  signal  is  then 
collected  by  a  boxcar  signal  averager  and  stored  in  a  computer  for  later 
analysis. 

1.4.2  Results 


1.4. 2.1  Vibrational  Relaxation.  A  three-photon  excitation  spectrum  of  hydrogen 
is  obtained  by  scanning  the  laser  and  measuring  the  total  H2(B)  fluorescence 
signal.  A  typical  spectrum  is  given  in  Fig.  19.  Alternatively  one  can  set  the 
laser  wavelength  to  excitate  a  single  rotational  transition  and  then  scan  the 
monochromator  to  obtain  a  resolved-fluorescence  spectrum.  The  spectrum  in 
Fig.  20  was  obtained  by  exciting  H2<B, v' =3 , J=2)  and  resolving  the  fluorescence 
signal  at  a  resolution  of  40A.  A  pressure  of  102  Torr  of  hydrogen  was  used. 

Our  computer  spectral  fitting  codes  were  then  used  to  generate  a  fit  to  the 
experimental  spectrum,  which  is  also  shown  in  Fig.  20.  All  features  in  the 
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Figure  19.  Three-photon  H2  excitation  spectrum.  H2(B,v'=3)  -  H2(X,v=0) 
80  Torr  of  82- 


WAVELENGTH  (nm) 

k-7743 

Figure  20.  Excitation  of  v'=3,  J'=2  transition  H2(B)  emission  spectrum. 
Spectrum  taken  with  40  A  resolution  and  102  Torr  of  H2. 


spectrum  could  be  assigned  to  v'=3,  the  vibrational  state  that  was  excited  by 
the  laser.  No  evidence  for  vibrational  relaxation  was  found.  Therefore 
vibrational  relaxation  is  not  an  important  product  channel  for  H2(B)  quenching. 


1.4. 2. 2  Quenching  Measurements.  To  determine  quenching  of  H2(B)  by  H2  the 
intensity  of  a  single  rotational  transition  was  measured  as  a  function  of  H2 
pressure  in  the  cell.  The  intensity  was  then  determined  by  integrating  the 
area  under  the  peak.  Use  of  peak  height  as  an  intensity  measure  was 
insufficient  for  two  reasons.  Pressure  broadening  of  the  peak  becomes 
appreciable  as  the  pressure  is  increased.  Therefore,  the  peak  height  does  not 
increase  as  rapidly  as  the  peak  area.  A  second  broadening  mechanism,  stark 
broadening  from  the  laser  field,  was  also  observed. 


1.4.3  Kinetic  Analysis 

The  production  and  depletion  of  H2(B)  can  be  summarized  in  three  equations: 

H2  +  h>j  (3  photons)  -*  H2(B,V'  ,J' )  k^  (19) 

H2(B,V',J')  -»  H2(X,V",J")  +  hv  (Lyman  Bands)  k2  (20) 

H2(B,V',J')  +  H2  -*  Products  (2H2,  H-Atoms,  H4*)  k3  .  (21) 

The  [ H2 ( B ) ]  is  in  steady  state  during  the  laser  pulse  since  the  radiative 
lifetime  is  ten  times  shorter  than  the  laser  pulse  duration: 


d[H2(B) ] 


^ -  =  k1[H2]  -  (k2  +  k3[H2J)  [HZ(B)1  =  0  .  (22) 


Substituting  [ H2 ( B ) J  =  I*/k2  and  rearranging,  one  obtains  the  equation 


[H-]  ,  k, 

-4“  =  n-  +  n~4”  lH2l  *  (23) 

I  kl  kl  k2  L 
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To  obtain  IC3  one  plots  [ H2 ] /I *  versus  [ H2 ]  where  the  slope/intercept  =  k.3/k.2- 
A  typical  plot  is  shown  in  Fig.  21.  This  plot  is  for  H2(B,v'=3,  J'=3) 
quenching  and  gives  a  quenching  rate  coefficient  of  13  ±  3  x  10~!0  cm^ 
molecule-*  s"!.  A  summary  of  all  the  quenching  rate  coefficients  determined  is 
given  in  Table  4.  There  does  not  appear  to  be  any  significant  difference  in 
quenching  coefficients  for  the  two  vibrational  levels  for  the  same  J  level. 
There  does  appear  to  be  some  difference  between  different  J  levels  but,  because 
of  the  uncertainty  in  the  measurements,  it  is  difficult  to  determine  if  a  true 
rotational  dependence  of  the  quenching  rate  coefficient  exists.  More 
measurements  are  needed  for  different  vibrational  levels  to  determine  if  there 
is  a  vibrational  dependence  of  the  quenching. 


[HU  (101 8  MOLECULES  cm'3) 


Figure  21.  H2(B,v'=3,  J'=3)  quenching  (k  =  13  ±  3  x  10"!*!  cm^  molecules"!  s"!) 


TABLE  4.  H2<B)  quenching  rate  coefficients. 


B 

D 

Transition 

k3(10-^  cm^  molecule--*  s-*) 

B 

0 

P(l) 

8.4  ±  1.8 

B 

m 

P(2) 

13  ±  3 

B 

R(l) 

15  +  6 

R(2) 

13  ±  3 

B 

R(3) 

5.5  +  1.5 

B 

B 

P  ( 2 ) 

6.6  +  1.4 

B 

R(l) 

13  ±  4 

■ 

D 

R(3) 

4.5  ±  1.0 
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2.  HIGH  SPIN  NITROGEN 

2.1  INTRODUCTION 

Experimental  investigations  of  "active"  nitrogen  have  long  postulated  the 
existence  of  a  high-energy  metastable  electronic  state  of  N£  responsible  for 
the  excitation  of  many  species.  Various  reports  indicate  the  state  has  a  low 
heterogeneous  deactivation  probability  on  glass  and  metal  surfaces  and  lives 
for  several  milliseconds  in  the  gas  phase.  Matrix  isolation  studies  have  shown 
that  chemiluminescence  from  "active"  nitrogen  trapped  at  10  to  35  K  may  radiate 
for  many  seconds. 25  The  existence  of  such  a  long-lived  metastable  state  would 
have  important  consequences  for  our  understanding  of  upper  atmospheric 
chemistry,  chemical  lasers,  energy  storage  and  high  voltage  switching. 

Recently,  several  theoretical  efforts  have  focused  on  the  stability  of  the 
"high  spin"  and  ^Eg+  electronic  states  of  N2.  These  states  arise  from  the 

pairing  of  two  NC^S)  atoms  to  form  an  N2  molecule  in  which  all  or  all  but  two 
of  the  electron  spins  remain  unpaired  in  the  resulting  molecular  orbitals. 

These  states  are  characterized  by  the  formation  of  a  weakly-bound  (several 
cm-l)  van  Der  Vaals  minimum  at  long  re  and  a  somewhat  more  strongly  covalently 
bound  minimum  at  shorter  re.  Since  these  configurations  are  the  lowest  energy 
of  each  of  the  spin  multiplicities,  they  should  be  stable  with  respect  to  decay 
via  dipole  radiation  and  perhaps  metastable  with  respect  to  collisional  deacti¬ 
vation.  Curve-crossing  to  the  state  has  been  postulated  by  Berkowitz, 
Chupka,  and  Kistiakowsky^  as  the  mechanism  for  predissociation  of  the  B^IIg 
state. 

Early  calculations  on  the  stability  of  the  ^E  state  were  conducted  by 
Vanderslice,  Mason,  and  Lippincott^  and  were  refined  by  Meador. 28  They 
obtained  13.4  and  14.7  cm-^-  respectively  for  the  depth  of  the  covalent  well. 
Subsequent  calculations  by  Ferrante  and  Stwalley29  and  Capitelli  and 
coworkers^O  obtained  well  depths  of  43  and  188  cm-!  respectively.  Recently, 
Partridge,  Langhoff,  and  Bauschlicher^l  performed  SDCI  ab  initio  calculations 
and  obtained  a  well  depth  of  21  cm-*. 
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The  potential  of  the  "’Eg+  state  was  initially  investigated  by  Mulliken,32  vho 
predicted  a  covalent  well  depth  of  -2300  cm-!.  Subsequent  calculations  by 
Krauss  and  Neumann^  predicted  a  shallower  potential  with  a  well  deptn  much 
closer  to  850  cm-^-.  Siegbahn^  calculated  that  De  was  approximately  1450  cm~l. 
Experimental  evidence  for  the  existence  of  a  stable  state  is  limited  to  the 
effects  it  may  have  on  the  predissociation  of  the  B-^Ilg  and 
proposition  advanced  by  Berkowitz  et  al.,26  that  the  ^E  state  is  responsible 
for  these  predissociations,  is  bolstered  by  the  evidence  of  Carroll, 35  who 
derived  a  well  depth  of  850  cm_3  from  predissociation  and  airglow  data. 

Recently,  Partridge  and  coworkers-^  performed  CASSCF/MRCI  ab  initio 

calculations  on  the  and  states  of  N2.  Their  calculations  for  the  ^E 

state  (designated  A')  revealed  a  covalent  well  depth  of  approximately  3450  cm_l 

and  a  van  Der  Uaals  minimum  with  a  binding  energy  of  47  cm~3.  A  barrier  of 

approximately  500  cm~l  with  respect  to  dissociation  to  form  two  N(^S)  atoms 

further  enhances  the  stability  of  the  covalently-bound  state.  Recombining 

N(^S)  atoms,  a  primary  mechanism  for  formation  of  the  A'-state,  must  traverse 

this  barrier  to  enter  the  covalently  bonded  well.  The  potential  supports  7 

bound  vibrational  levels  with  a  spacing  of  approximately  630  cm~3.  The 

potential  for  the  ^nu  state  (designated  C")  was  also  found  to  be  more  deeply 

bound  than  previously  estimated.  The  C"-state  well  depth  is  calculated  to  be 

10,505  cm-*  with  a  vibrational  spacing  of  approximately  925  cm~l.  A  potential 
% 

energy  diagram  of  these  states  is  shown  schematically  in  Fig.  22. 

Perhaps  the  most  important  result  from  the  Partridge  study  was  the  assignment 
of  the  previously  unassigned  Herman  Infrared  (HIR)  band  system  of  N2  to  the 
C"5nu"A'5Eg+  transition.  This  assignment  was  subsequently  confirmed  by  FTIR 
studies  of  Huber  and  Vervloet.3?  The  HIR  bands  are  observed  in  the 
650  to  900  nm  region  as  a  result  of  an  energy  pooling  reaction  between  two 
N2(A^Eu+)  molecules . 38-40  Emission  from  the  C^nu  and  B^IIg  states  is  also 
observed  from  the  energy  pooling  process.  The  calculated  radiative  lifetime  of 
the  C"-state  is  4.3  ps,  similar  to  the  B^ng-A^EU+  First-Positive  bands. ^  The 
calculated  lifetime  is  in  conflict  with  experimental  observations  that  the 
ratio  of  C"/C  state  emission  intensities  is  unchanged  over  a  pressure  range  of 


a^-IIg  states.  The 
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Figure  22.  Potential  energy  diagram  for  the  excited  quintet  and  triplet  states 
of  N£.  The  vibrational  levels  of  the  and  A'^Eg+  states  are 

indicated. 

a  few  Torr.  This  would  place  the  lifetime  of  the  C"-state  nearer  the  37  ns 
lifetime  of  the  C-state  and  raises  the  possibility  that  the  C"-state  derives 
its  lifetime  from  collisional  coupling  to  the  C-state. 

The  assignment  of  the  N2  HIR  bands  to  the  C"^ Hy-A' ^ £g+  transition  coupled  with 
the  3450  cm-1-  well  depth  establishes  a  means  of  production  and  indirect  evi¬ 
dence  for  the  stability  of  the  A'-state.  The  9.45  eV  electronic  term  energy  of 
the  A'-state  is  more  than  sufficient  to  excite  many  of  the  states  observed  in 
"active"  nitrogen.  However,  the  magnitude  of  the  intersystem  collisional 
transfer  processes  and  non-radiative  spin-orbit  and  spin-spin  couplings  to 
other  electronic  states  in  the  N2  manifold,  which  decrease  the  A'-state  life¬ 
time,  have  not  been  addressed.  The  experiments  reported  in  this  section  of  the 
report  employ  laser  induced  fluorescence  and  optogalvanic  techniques  to 
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quantitatively  explore  the  stability  of  the  A' -state  in  the  gas  phase.  Several 
methods  of  producing  the  A' -state  were  explored  as  part  of  this  effort.  The 
results  of  these  studies  are  described  in  the  sub-sections  which  follow. 

2.2  EXPERIMENTAL  APPROACH 


Ue  chose  to  try  three  different  techniques  to  form  the  N2  (A')  state.  The 
first  approach  was  a  DC  discharge  source  with  optogalvanic  detection.  Second 
was  an  RF  discharge  source  with  LIF  detection;  finally  a  DC  discharge  source 
with  LIF  detection.  Each  of  these  techniques  and  the  results  obtained  are 
described  in  the  following  three  sections. 

2.2.1  Optogalvanic  Spectroscopy 

In  optogalvanic  spectroscopy  a  tunable  laser  is  used  to  probe  the  spectral 
characteristics  of  atomic  or  molecular  species  generated  in  an  electric  dis¬ 
charge.  42*43  Two  conditions  are  required  to  generate  an  optogalvanic  signal. 
First,  the  laser  wavelength  must  equal  a  transition  of  a  species  present  in  the 
discharge.  Second,  the  absorption  of  energy  in  the  laser  beam  must  alter  the 
impedance  of  the  discharge.  There  is  general  agreement  that  almost  all  allowed 
transitions  will  cause  an  impedance  change  and  therefore  can  be  detected  with 
optogalvanic  spectroscopy.  There  is  less  agreement  as  to  the  cause  of  the 
impedance  change.  It  has  been  proposed  that  absorption  could  alter  the  concen¬ 
trations  of  rare  species  in  the  discharge  that  play  a  significant  role  in  the 
dynamics  of  the  discharge.  For  example,  absorption  could  result  in  the  pro¬ 
duction  of  easily  ionizable  states.  In  flames  it  has  been  observed  that  the 
ionization  potential  of  the  excited  state  does  influence  the  magnitude  of  the 
impedance  change.  However,  this  effect  has  not  been  observed  in  electric 
discharges.  Others  attribute  the  impedance  change  to  an  increase  in  the 
electron  energy.  This  could  cause  large  local  temperature  changes  to  which  the 
discharge  is  highly  sensitive. 

The  laser  induced  impedance  change  can  be  either  positive  or  negative.  The 
signal  is  measured  by  capacitively  decoupling  the  discharge  current  from  the 
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discharge  voltage.  The  resolution  of  the  optogalvanic  spectrum  is  determined 
by  the  laser  resolution.  The  detection  sensitivity  is  limited  by  the  electri¬ 
cal  noise  of  the  discharge.  For  atoms  a  detection  limit  of  10^  atoms  cm-^  is 
commonly  obtained.  The  sensitivity  for  molecules  tends  to  be  two  to  three 
orders  of  magnitude  lower  due  to  reduced  oscillator  strengths  and  distribution 
of  the  population  over  many  rotational  states. 

The  optogalvanic  spectrum  for  a  pure  N2  discharge  has  been  observed  between  570 
and  610  nm.^  The  spectrum  obtained  by  Feldmann  is  shown  in  Fig.  23.  Transi¬ 
tions  in  the  first  positive  (B^Ilg-A^£u+)  system  and  the  Ledbetter  Rydberg 
series  ( C4I  Ily-a"^  £gf)  were  observed.  Both  positive  and  negative  signals  were 
measured.  The  sign  of  some  of  the  transitions  changed  depending  on  the 
discharge  conditions. 
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Figure  23. 


Optogalvanic  spectrum  of  a  discharge  in  N2.  The  prominent 
bandheads  of  the  first  positive  system  are  indicated  by  the 
respective  vibrational  quantum  numbers.  (Ref.  44) 
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Experimental  Apparatus 


The  experimental  apparatus  is  shown  schematically  in  Fig.  24.  The  flow  reactor 
consists  of  a  1/2  in.  diameter  glass  tube  with  a  quartz  window  on  each  end. 

Gas  enters  through  a  port  at  one  end  of  the  tube  and  is  pumped  out  at  the  other 
end  by  a  27  cfm  mechanical  pump.  The  electrodes  are  made  of  tantalum  foil. 
Power  is  supplied  by  a  high  voltage  power  supply  in  series  with  a  current 
limiting  resistor.  The  discharge  typically  operates  at  500V  and  1  mA. 

Laser  excitation  is  provided  by  a  Quantel  system  consisting  of  a  YG580  series 
actively  Q-svitched  Nd:YAG  laser  pumping  a  Quantel  dye  laser.  The  laser  dye 
LDS698  is  used  to  obtain  the  needed  excitation  wavelength.  Typical  output 
energies  of  30  mJ  per  pulse  are  obtained  from  this  laser  of  a  rate  of  10  pulses 
per  second.  The  optogalvanic  signal  is  collected  using  a  gated  boxcar  integra¬ 
tor.  The  averaged  signal  is  then  stored  on  a  chart  recorder  or  computer.  A 
0.35m  McPherson  monochromator  with  a  cooled  R943-02  Hamamatsu  photomultiplier 
tube  detector  is  used  to  obtain  a  spectrum  of  the  discharge  emission.  In  this 
way  we  can  verify  that  N2<A)  is  being  produced  and  energy  pooling  to  form  the 
HIR  bands  is  occurring. 
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Results 


Discharge  Emission  Spectrum 

A  portion  of  the  discharge  emission  spectrum  recorded  with  the  monochromator  is 
shown  in  Fig.  25.  Observation  of  the  (3,1)  and  (2,0)  bands  of  the  HIR  system 
is  clear  evidence  that  N2(A)  energy  pooling  is  occurring  in  the  discharge. 

Also  observed  is  emission  from  the  first  positive  system,  ^(B^ng-A^Ey*)  and 
argon  metastables.  The  two  HIR  bands  shown  have  no  overlap  with  the  first 
positive  system  and  therefore  can  be  observed  cleanly.  Many  of  the  other  HIR 
bands  are  obscured  under  more  intense  first  positive  emission.  For  this  reason 
the  region  of  the  (3,1)  and  (2,0)  bands  was  chosen  as  the  region  to  probe  for 
optogalvanic  signal  from  the  laser. 

The  rate  of  C"-state  formation  may  be  determined  directly  via  absolute 
intensity  measurements  of  the  HIR  bands.  The  absolute  intensities  of  spectra 
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Figure  25.  Spectrum  of  the  nitrogen  Herman  infrared  and  first-positive  system 
systems  excited  in  the  energy  pooling  of  N2(A)  in  a  nitrogen/argon 
discharge. 
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such  as  Fig.  25  can  be  integrated  to  obtain  the  total  intensity  of  N2(C"-A'), 

2-0  and  3-1  bands.  Other  vibrational  levels  of  the  C"  state  are  not  produced 
from  energy  pooling.  Since  all  emission  for  the  C"  state  does  not  occur  from 
these  transitions,  the  2-0  and  3-1  band  intensities  must  be  corrected  for  the 
fluorescence  branching  to  other  lover  state  vibrational  levels  in  order  to 
obtain  the  total  emission  rates  from  these  levels.  Finally,  the  total  emission 
intensities  from  these  levels  must  be  summed  to  get  the  total  HIR  formation 
rate.  The  radiative  transition  probabilities  from  Langhoff  and  coworkers^  are 
employed  to  calculate  the  fluorescence  branching  ratios. 

The  N2(C")  formation  rate  based  on  the  intensities  of  the  2-0  and  3-1  HIR  bands 
is  7.0  x  loll  molecules  cm^  s~l.  This  rate  is  also  a  lover  limit  on  the  rate  of 

N2(A')  formation,  since  the  A'  state  may  be  produced  by  other  mechanisms  besides 

radiative  cascade  from  the  C"  state.  The  N2(A')  population  may  be  determined  if 
information  about  the  A'  state  loss  rates  are  available.  Such  processes  as 
collisional  quenching  and  collision-induced  predissociation,  curve  crossing, 
radiation,  and  collisional  quenching  on  the  vails  must  be  considered.  An  upper 
limit  to  the  population  may  be  obtained  by  considering  only  vail  losses. 

The  time  required  for  a  molecule  to  diffuse  from  the  center  of  a  cylinder  to  the 

vails  is  given  by 


(24) 


where  r  is  the  flow  tube  radius  and  D  is  the  diffusion  coefficient  ( —30  cm^s-* 
at  6  Torr).  The  flov  tube  radius  is  1.27  hence  the  characteristic  diffusion 
time  is  22.4  ms.  This  is  comparable  to  the  residence  time  of  the  species  in 
the  flov  reactor  and  hence  a  steady  state  solution  can  be  used  to  obtain  the 
N2(A')  number  density: 

d[N  (A')1 

dt  -  °  =  krF2(C,,)]  -  kPF2(A,)]  •  (25) 
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In  this  expression,  kr(N2(C")|  is  the  radiative  rate  of  ^(C"),  Ijjjr,  and  kp  is 
the  inverse  diffusion  time.  Hence,  Eq.  (25)  can  be  rewritten 


P2<*' >]  -  ^  '  Wb 


(26) 


Solving  Eq.  (26)  using  the  HIR  intensity  and  N2(A')  diffusion  time  quoted  above 
gives  an  N2(A' )  density  of  approximately  1.6  x  10*0  molecule  cm-^. 

Optogalvanic  Spectrum 

The  excitation  laser  was  scanned  over  the  range  695  to  716  nm  and  the  resulting 
optogalvanic  spectrum  was  recorded.  Many  transitions  were  observed.  A  second 
optogalvanic  spectrum  of  a  low  pressure  argon  discharge  was  recorded  to  aid  in 
distinguishing  between  argon  and  nitrogen  transitions.  Almost  all  of  the 
transitions  observed  in  the  ^/Ar  discharge  are  also  present  in  the  pure  Ar 
discharge.  A  portion  of  the  spectrum  for  each  discharge  is  shown  in  Figs.  26 
and  27.  The  few  remaining  unassigned  lines  cannot  be  readily  assigned  to  any 
nitrogen  transition.  They  could  belong  to  argon  transitions  that  are  too  weak 
under  the  low  discharge  conditions  of  the  pure  Ar  spectrum. 

It  was  impossible  to  record  a  pure  Ar  spectrum  at  the  same  pressure  as  the 
N2/Ar  discharge  because  the  discharge  noise  was  as  large  as  the  signal  at 
higher  pressures  of  pure  argon.  The  discharge  noise  for  the  ^/Ar  discharge 
actually  got  smaller  when  the  pressure  was  increased  to  above  4  Torr.  Spectra 
were  recorded  under  different  discharge  conditions  to  determine  the  effect  on 
the  peak  intensities.  Three  spectra  taken  under  different  conditions  are  shown 
in  Figs.  28  to  30.  The  relative  intensities  of  the  three  peaks  are  dramati¬ 
cally  different  for  the  three  different  discharge  conditions.  The  third  peak 
appears  as  a  doublet  in  the  pure  Ar  spectrum  but  as  only  a  single  peak  in  the 
N2/Ar  discharge  spectra.  Considering  the  intensity  fluctuations  that  are 
observed  as  a  function  of  discharge  condition  it  is  possible  that  all  the 
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Figure  28.  Optogalvanic  spectrum  of  a  nitrogen/argon  DC  discharge  at  4.50  Torr. 
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Figure  29.  Optogalvanic  spectrum  of  a  nitrogen/argon  DC  discharge  at  2.30  Torr. 
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Figure  30.  Optogalvanic  spectrum  of  an  argon  DC  discharge  at  0.97  Torr. 

unassigned  lines  are  argon  transitions.  There  is  no  real  evidence  to  suggest 
that  any  of  the  transitions  observed  belong  to  N2. 

Discussion 


The  observed  emission  spectrum  confirms  that  our  discharge  was  operating 
properly.  The  observation  of  HIR  emission  indicates  that  the  A'  state  was 
present  in  the  discharge.  Our  calculations  indicate  the  A'  state  may  be 
present  at  densities  approaching  1.6  x  10^  molecule  cm"3.  However  we  were 
unable  to  observe  any  optogalvanic  signal  from  any  excited  N£.  There  are 
several  reasons  that  this  might  occur.  Our  discharge  source  is  homemade  and 
was  much  noisier  than  expected.  Achieving  a  low  noise  discharge  for  optogal¬ 
vanic  spectroscopy  appears  to  be  inconsistent  with  optimal  production  of  the 
HIR  band  system.  Every  gas  discharge  has  a  different  noise  frequency  spectrum 
and  methods  for  reducing  noise  also  vary  for  each  type  of  discharge. 
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Previous  reports  on  the  sensitivity  of  optogalvanic  spectroscopy  with  ^eta- 
stable  N2(A^Eu+)  ind  icate  that  densities  as  low  as  10®  to  10^  molecules  cm-^ 
may  be  observed. ^4  Uhen  we  attempted  to  excite  these  transitions  in  our 
system,  we  were  unable  to  observe  impedence  changes  corresponding  to  the 
N2(B^ng  -  A^EU+)  First  Positive  Bands.  We  estimate  the  ^(A^E)  density  to  be 
approximately  3  x  10^  molecule  cm-^  in  this  system.  Hence,  a  better  estimate 
of  the  detection  sensitivity  under  our  flow  conditions  is  approximately 
10*0  molecule  cm~3.  The  calculations  of  Partridge  et  al.36  show  that  B-A  and 
C"-A'  have  similar  oscillator  strengths.  Our  calculations  indicate  that  even 
modest  loss  processes  would  reduce  the  N2(A' )  number  densities  expected  in  our 
system  below  the  rough  detection  limit  for  these  particular  experiments. 

Hence,  it  is  quite  likely  that  the  A'-state  is  produced  quite  efficiently  under 
our  conditions  but  that  optogalvanic  spectroscopy  does  not  have  sufficient 
sensitivity  to  detect  it  at  the  levels  present  in  the  flow  reactor. 

2.2.2  Laser-Induced-Fluorescence  RF  Discharge  Source 

Another  possible  method  to  produce  the  A'-state  is  by  N-atom  recombination. 

The  calculations  of  Partridge  et  al.^l  show  that  recombination  of  N-atoms  on 
the  ^Z  potential  is  enhanced  at  lower  temperatures.  We  produced  N-atoms  with 
an  RF  discharge  source  in  a  supersonic  nozzle.  The  cooling  expected  during  the 
expansion  was  expected  to  enhance  the  A'-state  formation.  N-atom  recombination 
occurs  at  the  throat  of  the  nozzle.  We  used  laser-induced-fluorescence  (LIF) 
on  the  HIR  bands  to  probe  for  the  presence  of  N2(A')  in  the  gas  beam  downstream 
of  the  nozzle.  The  following  sections  contain  a  description  of  the  RF 
discharge  source  and  the  modifications  necessary  to  incorporate  the  source  in 
our  existing  vacuum  chamber. 

The  RF  discharge  source  was  based  on  the  design  of  Lee  and  coworkers. ^  He 
developed  a  high  pressure,  RF  discharge  nozzle  beam  source  capable  of  producing 
a  very  intense  (>  10*®  atoms  sr-*  s-*)  supersonic  beam  of  oxygen  atoms.  Two 
important  developments  lead  to  the  achievement  of  such  an  intense  source. 

First,  an  efficient  impedance  matching  scheme  was  devised  for  coupling  the  RF 
power  to  the  oxygen-rare  gas  mixtures.  Second,  techniques  were  developed  for 
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localizing  the  discharge  directly  behind  the  orifice  of  the  nozzle.  The  second 
development  is  quite  important  since  atoms  will  recombine  quite  efficiently 
inthe  slow  gas  flow  downstream  of  the  actual  discharge.  Extension  of  the 
discharge  to  the  nozzle  throat  prevents  significant  recombination  from 
occuring  prior  to  the  expansion. 

A  cross-sectional  view  of  the  source,  with  an  enlarged  view  of  the  nozzle  tip, 
appears  in  Fig.  31.  The  nozzle  tube  is  made  of  quartz.  Lee  tested  two 
materials,  quartz  and  alumina.  These  were  chosen  because  of  their  relatively 
high  melting  points.  Alumina  tubes,  however,  were  found  to  be  unsatisfac tory . 
They  tended  to  crack  when  subjected  to  large  thermal  gradients  and  arcing 
occurred  between  the  outer  surface  of  the  tube  and  the  RF  coil  whenever  the  two 
came  into  physical  contact.  Therefore  quartz  was  chosen  as  the  nozzle  material 
despite  being  more  difficult  to  fabricate  into  a  nozzle. 

In  Lee's  final  design  it  was  found  necessary  to  incorporate  a  water  cooling 
jacket  to  prevent  the  quartz  orifice  at  the  tip  from  melting  and  enlarging 
during  the  expansion  of  the  hot  gas  mixture.  Quartz  is  a  poor  thermal  conductor 
so  it  was  important  that  the  water  coolant  flow  be  as  close  to  the  orifice  as 
possible.  With  the  water  jacket  the  orifice  was  found  to  be  stable  with  respect 
to  enlargement.  It  is  necessary  to  use  low  electrical  conductivity  water  as 
coolant  to  reduce  excessive  RF  power  loss  to  the  water.  Lee  found  an  added 
benefit  of  the  jacket  was  the  elimination  of  arcing  between  the  plasma  tube  and 
the  RF  coils. 

In  our  system  the  RF  power  is  supplied  by  an  ENI  Power  Systems  ACG5  plasma 
generator  which  can  deliver  a  maximum  output  of  500V..  Impedance  matching  is 
accomplished  by  an  Nye  Viking  MB-IV-A  impedance  matching  network.  Impedance 
matching  is  necessary  since  the  effective  plasma  impedance,  which  varies  as  a 
function  of  plasma  condition,  is  much  larger  than  the  508  output  impedance  of 
the  electronics.  A  schematic  of  the  RF  discharge  circuit  is  given  in  Fig.  32. 
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a.  Variable  ground  tap 

b.  Svagelock  reducer 

c.  Water-cooled  quartz  nozzle 

d.  Coupling  coil 

e.  Quartz  support  rod  for  water'  inlet/output  assembly 

f.  Svagelock  union  joining  copper  and  polyflo  tubing 

g.  Cajon  VCO  fitting 

h.  Kovarpyrex  section 

i.  Graded  seal,  Pyrex  to  Vycor 

j.  Aluminum  support  block  to  electrical  ground 

k.  Boron  nitride  skimmer,  0.88  mm 

l.  4x  enlargement  of  nozzle  tip 


Figure  31.  Sectioned  view  of  the  internal  source  components  (arrows 

indicate  low  conductivity  water  flow  direction.  (Ref.  46)) 
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Figure  32.  Schematic  of  radio  frequency  source. 

A  variable  ground  tap  on  the  tank  coil  is  used  to  produce  a  large  stepdovn  of 
the  plasma  impedance,  as  seen  by  the  RF  electronics,  ideally  to  50Q.  With  this 
RF  circuit  it  is  possible  to  operate  routinely  at  a  standing  wave  ratio  (SWR) 
of  less  than  1.05:1  for  02/rare  gas  mixtures.  For  N2  we  found  that  the  coupl¬ 
ing  was  less  efficient  but  that  a  SWR  of  less  than  2:1  should  be  attainable. 

To  localize  the  plasma  directly  behind  the  nozzle  orifice  an  electrically 
grounded  block  of  aluminum  (Fig.  31, J)  is  located  around  the  nozzle  tip.  At 
high  pressures  coupling  occurs  between  the  front,  small  diameter  tank  coil 
turns  and  the  aluminum  block.  As  this  coupling  becomes  stronger  the  plasma 
localizes  towards  the  front  of  the  discharge  tube.  The  tank  coil  (Fig.  31, D) 
is  differentially  wound  to  further  localize  the  plasma.  The  larger  turns 
decouple  from  the  plasma  at  higher  pressures  and  increase  the  energy  density  of 
the  localized  plasma.  It  is  important  to  note  that  the  entire  tank  circuit,  as 
shown  in  Fig.  32,  is  floated  with  respect  to  electrical  ground.  This  is 
required  if  the  RF  coupling  and  the  localization  are  to  be  effective. 
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Chamber  Modifications 


Our  chamber  used  in  the  H4*  experiments  consisted  of  a  10  in.  six-way  cross 
with  a  smaller  detection  chamber  attached  at  one  end.  Both  chambers  are  pumped 
by  a  single  diffusion  pump.  We  modified  the  system  to  create  two 
differentially  pumped  chambers.  A  stainless  steel  plate  with  a  small  hole  was 
inserted  between  the  two  chambers  and  a  second  diffusion  pump  was  added  to  the 
small  detection  chamber.  A  boron  nitride  skimmer  was  mounted  on  this  plate  on 
the  large  chamber  side. 

The  RF  discharge  source  is  mounted  on  a  10  in.  conflat  flange  which  mounts  on 
the  large  chamber  directly  opposite  the  detection  chamber.  Both  RF  and  water 
vacuum  feedthroughs  to  the  discharge  source  have  been  made.  Mounting  brackets 
to  hold  the  nozzle  in  place  after  alignment  have  also  been  made.  Wherever 
possible  non-metallic  materials  have  been  used  to  minimize  coupling  of  the  RF 
to  mounting  and  support  brackets. 

Since  we  were  using  LIF  detection,  it  was  necessary  to  assemble  Brewster  angle 
windows  with  light  baffles  for  the  laser  light  to  pass  through.  The 
photomultiplier  (PMT)  used  to  detect  the  LIF  signal  was  equipped  with  a  band¬ 
pass  filter  to  isolate  the  HIR  band  transition  of  interest.  However,  the  laser 
light  is  edge  of  the  filter  of  the  filter  transmission  curve  wavelength  of  the 
quite  close  to  the  so  scattered  light  must  be  carefully  eliminated.  A  1:1 
imaging  lens  system  was  used  to  collect  the  fluorescence  focus  it  on  the  PMT, 
and  spatially  define  the  collection  region  to  minimize  scattered  light. 

Unfortunately,  our  chamber  is  considerably  smaller  than  Lee's,  making  for  a 
very  close  fit.  The  proximity  of  the  RF  coils  to  the  chamber  walls  allows  the 
RF  radiation  to  couple  to  the  chamber  instead  of  the  discharge  tube.  Through¬ 
out  our  testing  of  the  RF  source  we  were  unable  to  sustain  a  discharge  at  the 
nozzle  tip.  We  therefore,  replaced  the  nozzle  assembly.  The  new  nozzle  had 
the  RF  discharge  just  outside  the  vacuum  chamber.  The  nozzle  extended  into  the 
chamber  several  inches  but  remained  out  of  the  detector  f ield-of-view. 
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Laser-induced- fluorescence  was  used  to  probe  for  both  HTR  and  the  first  positive 
bands  of  N2.  No  signal  vas  observed.  The  same  experiment  was  also  performed 
with  the  RF  discharge  replaced  with  a  microwave  cavity.  Again  no  positive 
results  were  obtained.  It  was  then  decided  to  concentrate  the  remainder  of  the 
program  on  the  DC  discharge  flow  tube  source  where  we  were  having  some  success. 

2.2.3  Dishcarge  Flow  Reactor/LIF  Studies 

Our  failure  to  observe  the  presence  of  any  ^(A')  using  the  optogalvanic 
detection  system  or  the  RF  excited  source  prompted  us  to  re-examine  the  use  of 
LIF  with  a  discharge  flow  reactor.  Originally,  we  felt  that  the  N2(A') 
molecule  might  be  efficiently  quenched  in  the  gas  phase.  As  a  result  of  this 
assumption,  we  had  sought  to  detect  the  A'-state  within  the  discharge  region 
(using  optogalvanic  spectroscopy)  and  in  a  collision-free  molecular  beam.  The 
use  of  the  discharge  flow  reactor  allows  us  to  better  quantify  the  A'-state 
production  role  through  observation  of  the  HIR  band  emission.  We  are  also  able 
to  produce  and  detect  LIF  over  a  larger  collection  volume  in  the  flow  reactor 
than  in  a  molecular  beam.  Finally,  simultaneous  emission  and  LIF  detection  was 
possible  in  the  discharge  flow  reactor.  The  results  of  this  final  set  of 
experiments  are  outlined  below. 

Our  approach  to  determining  the  collisional  stability  of  the  A'^Sg+  state  o  f  N2 
was  as  follows.  We  employed  energy  pooling  of  two  ^(A^Ey*)  molecules  to  pro¬ 
duce  the  C^llg  state.  Radiative  cascade  through  the  C"-A'  HIR  bands  results  in 
production  of  the  A'-state.  The  production  rate  can  be  determined  absolutely 
via  quantitative  observation  of  the  HIR  band  emission  intensity.  Laser-induced 
fluorescence  (LIF)  on  the  C"-A'  transition  is  employed  to  probe  for  the  presence 
of  the  A'-state.  Calibration  of  the  LIF  detection  system  is  accomplished 
approximately  by  comparison  with  LIF  signals  from  the  nearby  B^ng-A^EU+  transi¬ 
tion.  The  N2(A^E)  population  may  be  inferred  from  HIR  band  emission  intensities 
using  the  HIR  energy  pooling  rate  coefficient  or  directly  using  N2(A-X) 
Vegard-Kaplan  band  emission  intensities.  A  simple  kinetic  model,  tied  to 
absolute  population  measurements,  is  used  to  predict  the  A'-state  population. 

The  HIR  band  LIF  intensities  were  compared  with  intensities  predicted  by  the 
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calibrated  computer  model.  The  gas  phase  ^(A')  lifetime  required  to  explain 
the  observed  signal  level  vas  then  determined  from  the  model.  The  experiment  is 
described  in  more  detail  below. 

The  discharge  flow  apparatus  is  shown  schematically  in  Pig.  33.  The  apparatus 
consists  of  an  ^(A-^I)  generator,  a  short  flow  section,  diagnostics  section, 
and  pumping  facility.  The  flow  section  vas  provided  to  allow  the  N2<A')  to 
build  to  maximum  concentration,  to  reduce  scattered  light  from  the  discharge, 
and  to  add  various  reagents.  The  diagnostics  section  provided  for  observation 
of  passive  emission  and  LIF  at  the  same  point  in  the  flow  reactor.  Passive 
emission  was  detected  using  a  0.35m  monochromator  (McPherson  270)  equipped  with 
a  cooled  GaAs  PMT  (HTV  R943-02)  operated  in  photon  counting  mode.  A  rate  meter 
on  the  photon  counter  (Pacific  Photometries)  was  interfaced  to  a  laboratory 
computer  for  signal  acquisition.  The  absolute  sensitivity  of  the  detection 


Figure  33.  Schematic  diagram  of  the  discharge  flow  reactor  and  LIF  detection 
system  used  in  the  experiments. 
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system  was  calibrated  using  the  well  known  A/NO  air  afterglow  in  the  400  to 
900  nm  region. 40,46  The  absolute  calibration  was  extended  to  200  nm  by 
comparison  with  deuterium  and  quartz  halogen  irradiation  standards. 

A  Nd:YAG  pumped  tunable  dye  laser  (Quantel  YG581C/TDL50)  operating  with  LDS  698 
dye  was  used  to  excite  fluorescence  from  the  N2  A' -state  and  A-state  in  the 
diagnostics  section.  The  section  is  internally  blackened  with  a  Teflon  coating 
and  equipped  with  long  baffle  arms  to  reduce  scattered  light.  The  LIF  is  detec¬ 
ted  with  a  baffled  and  filtered  PMT  (HTV  R955).  An  f/2  lens  was  used  to  image 
fluorescence  on  to  the  PMT  photocathode.  The  combination  of  a  Schott  RG-9  color 
filter  and  750  nm  long-pass  interference  filter,  along  with  the  PMT  long-wavele¬ 
ngth  cut-off,  limited  the  wavelength  range  of  the  detection  system  to  750  to 
900  nm.  Approximately  50  percent  of  the  emission  from  N2(B,  v=3)  and  N2(C", 
v=2),  the  levels  pumped  in  the  LIF  detection  experiments,  occurs  in  this  wavel¬ 
ength  region.  Signals  from  the  PMT  were  fed  to  a  high-speed  current  amplifier 
(Ithaco  1211)  and  then  to  a  boxcar  signal  averager.  The  LIF  signals  were  recor¬ 
ded  with  a  laboratory  computer  and  corrected  for  variations  in  laser  intensity. 

The  N2(A^E)  used  in  this  study  was  produced  by  energy  transfer  from  metastable 
Ar(^Po>2)  to  n2*  This  technique  has  been  described  extensively  in  the  litera¬ 
ture.  47,48  The  metastable  Ar  was  produced  by  passing  Ar  though  a  low  power  DC 
discharge  (500  VDC,  2  mA)  at  a  flow  rate  of  approximately  4000  untole  s_*. 

Small  flows  of  N2  (10  to  250  ymole  s-*)  were  added  downstream  of  the  discharge. 
In  some  experiments  a  5  percent  Ar  in  He  mixture  was  passed  through  the 
discharge  in  order  to  use  He  as  the  buffer  gas.  A  small  flow  of  O2  could  be 
injected  downstream  of  the  ^(A^E)  generator  in  order  to  quench  the  N2(A). 

This  was  done  occasionally  in  order  to  eliminate  all  chemiluminescence  and 
check  for  the  presence  of  scattered  light  from  the  discharge.  In  all  cases 
flows  were  measured  using  calibrated  mass  flow  meters.  Cell  pressures  were 
measured  using  a  calibrated  capacitance  manometer.  All  reagent  gases  were  UHP 
(99.999)  and  were  used  without  further  purification.  The  flow  reactor  was 
pumped  by  two  17  cfm  and  one  27  cm  mechanical  vacuum  pumps  connected  in  a 
parallel  manifold.  In  this  configuration  flow  velocities  in  excess  of 
3700  cm  s-l  at  pressures  of  2  Torr  were  attainable. 
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2.2.4  Results 


Initial  experiments  were  conducted  using  Ar  as  a  buffer  gas  at  pressures  of 
approximately  4  Torr.  Observation  of  N2  HIR  emission  shoved  that  the  intensity 
was  maximized  under  these  conditions  at  N2  mole  fractions  of  approximately 
0.05.  An  emission  spectrum  of  the  region  from  610  to  820  nm,  recorded  at  a 
point  5  ms  downstream  of  the  N2<A)  generator,  is  shown  in  Fig.  34.  The 
spectrum  includes  the  data  and  a  fit  to  the  data  indicating  the  presence  of 
both  the  N2  (C"-A')  HIR  bands  and  the  N2(B-A)  First  Positive  bands.  The 
transition  probabilities  calculated  by  Partridge  et  al.^6  were  used  in 
calculating  the  fitted  spectrum.  It  is  evident  that  these  parameters  do  quite 
a  good  job  of  reproducing  the  experimental  spectrum.  The  spectrum  of  the  C"-A' 
transition  is  dominated  by  emission  from  vibrational  levels  2  and  3.  Previous 
studies  by  Piper^O  showed  that  these  levels  are  produced  nearly  exclusively  by 
the  energy  pooling  reactions  involving  combinations  of  N2(A,  v=0)  and 
N2(A,  v=l ) 


610  660  710  760  810 


WAVELENGTH  (nm)  M7SSa 

Figure  34.  Visible  emission  spectrum  of  N2  resulting  from  the  N2(A)  energy 
pooling  reaction.  The  important  C"-A'  and  B-A  transitions  are 
indicated . 
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A  lower  limit  on  the  N2(A')  production  rate  is  given  by  the  C"-A'  emission 
intensity.  This  emission  intensity  can  also  be  used  to  determine  the  N2(A) 
concentration  at  that  point  in  the  flow  reactor.  The  total  N2(A')  concentra¬ 
tion  at  any  point  in  the  flow  reactor  is  the  integral  of  N2(A')  production  and 
loss  processes  from  the  N2(A)  generator  to  that  location.  Initially,  we  assume 
that  only  quenching  at  the  walls  of  the  flow  reactor  is  responsible  for  the 
loss  of  N2(A')*  The  production  and  loss  of  N2(A)  and  N2(A')  can  be  described 


approximately  by  the  set  of  reactions: 

N2(A)  +  N2(a>  N2(C")  +  n2<x>  (27) 
N2(A)  +  N2(A)  ->N2(B,C)  +  N2(X)  (28) 
N2(A)  +  wall  -»  N2(X)  (29) 
N2(C")  -»N2(A')  +  hw(HIR)  (30) 
N2(A')  +  wall  ->  N2(X)  .  (31) 


Other  coll isional  and  non-collisional  pathways  for  deactivation  of  N2(A')  can 
be  described  by  inclusion  of  an  additional  pseudo  first-order  loss  term: 

N2(A')  -»  products  (32) 

The  radiative  rate  of  the  C"-state  is  sufficiently  large  so  that  the  C"-state 
concentration  may  be  considered  to  be  in  steady  state  in  the  field  of  view  of 
the  detection  system.  Hence,  the  rate  of  Eq.  (27)  can  be  set  equal  to  the  rate 
of  Eq.  (30): 

ki[N2(A)!2  =  k4[N2(C')l  =  IHIR  (33) 

where  IfjlR  is  simply  the  total  Herman  Infrared  Band  volume  emission  rate. 

Thus,  the  N2(A)  concentration  may  be  determined  from  the  HIR  band  absolute 
emission  rate: 

[N2(A)1  =  (I(HIR)/k1) 1/2  .  (34) 
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The  energy  pooling  rate  coefficients  for  the  formation  of  the  N2(C")  state  have 
been  measured  by  Piper. ^0  The  rate  coefficients  and  the  vibrational  levels  of 
the  C"-state  that  are  excited  depend  on  the  degree  of  vibrational  excitation 
present  in  the  N2<A).  Under  our  experimental  conditions  the  N2(A)  is  produced 
primarily  in  v=0.  Hence,  ve  have  used  Piper's  rate  coefficient  for  pooling  of 
two  N2<A,  v=0)  molecules,  8.1  +  2.3  x  10“^  cm^  molecule-!  s-!,  in  determining 
the  N2(A)  number  density.  Thus,  under  typical  experimental  conditions  the  HIR 
band  emission  measurement  gives  an  N2<A)  number  density  of  9  x  10^  molecules 
cm-3  5  ms  downstream  of  the  generator. 

This  estimate  of  the  N2(A)  number  density  is  confirmed  by  direct  observation  of 
the  N2<A-X)  Vegard-Kaplan  Bands.  Figure  35  shows  a  spectrum  of  the  region  from 
200  to  400  nm.  The  fit  to  the  spectrum  clearly  identifies  features  from 
N2<A,  v=0,1)  and  N2(C,  v=0-3).  The  N2<C)  emission  is  produced  from  the  N2(A) 
energy  pooling  reaction.  The  emission  from  the  NO(A-X)  transition  evident  in 


Figure  35.  UV  emission  spectrum  of  N2  from  N2(A)  energy  pooling.  Features  of 
the  A-X  and  C-B  transitions  are  indicated.  Also  present  are 
features  from  the  NO(A-X)  and  OH(A-X)  transitions  resulting  from 
the  presence  of  low-level  impurities. 
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the  spectrum  is  due  to  excitation  of  small  amounts  of  discharge-produced  NO  by 
N2<A).  Also  present,  but  not  identified  in  the  spectrum,  is  emission  from  the 
OH(A-X)  transition  which  is  also  due  to  a  low-level  impurity.  The  direct 
determination  of  N£(A)  concentrations  from  the  A-X  band  intensities  in  general 
indicated  N£(A)  concentrations  consistent  with  HIR  emission  intensity 
determinations. 

The  decay  of  N2(A)  in  the  flow  reactor  can  be  modeled  to  determine  the  initial 
N2(A)  concentration  at  the  exit  of  the  generator.  This  concentration  must  be 
determined  in  order  to  calculate  the  total  N2(A')  production  up  to  the 
observation  point.  The  decay  of  N2(A)  is  given  by  the  expression: 

d[N  (A) ] 

—ft -  =  (k1+k2)[N2(A)]Z  -  k3[N2(A)]  (35) 

where  the  factor  of  two  usually  found  in  the  rate  expression  reflecting  the 
loss  of  two  N2<A)  molecules  per  pooling  event  has  been  excluded.  This  is  done 
to  account  for  the  fact  that,  under  most  experimental  conditions,  one  N2(A) 
molecule  is  regenerated  from  the  pooling  reaction  via  radiative  cascade  through 
the  C3nu  and  B3IIg  states.  The  rate  coefficient  for  wall  loss  in  a  diffusion 
dominated  regime,  k3,  is  given  by:^^ 

k3  =  3.66  D/a2  (36) 

where  D  is  the  diffusion  coefficient  of  N2(A)  in  the  buffer  gas  and  a  is  the 
tube  radius.  The  diffusion  coefficient  can  be  calculated  approximately  from 
kinetic  theory3®  to  be  160  cm2  s-3  at  a  pressure  of  1  Torr  and  scales  inversely 
with  pressure.  Under  our  conditions  k3  typically  has  a  value  of  80  s-3. 

The  value  of  k2,  the  energy  pooling  rate  coefficient  for  production  of  all 
products  other  than  the  C"-state,  is  estimated  by  Piper^®  to  be  approximately 
3  x  10~3®  cm3  molecule-3  s~3.  Thus,  it  is  clear  that  at  these  N2(A)  concentra¬ 
tions  wall  loss  is  primarily  responsible  for  the  loss  of  N2(A).  Equation  (36) 
can  be  solved  analytically  by  neglecting  pooling  or  exactly  using  a  simple 
numerical  model  to  yield  the  ^(A)  concentration  at  the  generator  exit  port. 
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The  equation  governing  the  production  and  loss  of  N2(A')  is  given  by: 


d[Nz(A')] 


k4[N2<C")l  -  {k5  +  k6)[N2(A')l 


where  the  value  of  k5,  wall  loss  rate  coefficient  for  the  A'-state,  is  assumed 
to  be  the  same  as  for  the  N2(A)  state  and  k4[N2(C")]  is  shown  in  Eq.  (33)  to  be 
the  total  HIR  band  emission  intensity.  By  substitution  of  Eq.  (33)  into 
Eq.  (37)  one  can  relate  the  expected  A'-state  concentration  to  the  N2(A) 
concentration: 


d[N2(A')J 


=  k1[N2(A)]Z  -  {k5  +  k6}[N2(A')] 


Given  the  initial  N2(A)  concentration  determined  from  Eq.  (36)  the  concentra¬ 
tion  of  N2(A' )  can  be  calculated  at  any  point  in  the  flow  reactor  using  a 
simple  numerical  model.  It  was  initially  assumed  that  other  loss  processes 
were  negligible  (kg  =  0).  The  model  predicted  that  N2(A')  concentrations  on 
the  order  of  4  x  10^  molecules  cm-^  could  be  expected  in  the  observation 
region,  5  ms  downstream  of  the  injector. 

Approximate  calibration  of  the  LIF  detection  system  was  accomplished  by 
observing  the  LIF  produced  by  pumping  the  ^(B-^IIg-A^ Eu+,  3-0)  band.  The  origin 
of  the  band  is  at  685.83  nm.  This  is  a  spectral  region  where  no  absorption  is 
expected  from  the  ^(C'-A' )  transition.  A  typical  spectrum  of  LIF  from  the  3-0 
band  is  shown  in  Fig.  36.  The  transition  is  a  ^II(Case  a)-^E  transition  with  27 
possible  branches.  Several  of  the  main  branches  and  satellite  branches  are 
indicated  in  the  figure.  The  signal  intensity  at  the  Pjj  branch  head  was 
chosen  for  determining  the  detection  sensitivity.  The  absolute  ^(A^E,  v=0) 
concentration  was  determined  by  the  techniques  described  above. 

Several  issues  must  be  considered  in  using  this  method  to  determine  the 
sensitivity  of  the  LIF  technique  for  detecting  the  A'-state.  The  primary  issue 
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WAVELENGTH  (nm) 

Figure  36.  LIF  excitation  spectrum  of  the  B^ng-A^J^*  3,0  band  transition 

observed  in  the  discharge  flow  reactor.  The  and  branch 
band  heads  are  indicated  in  the  figure. 

is  the  actual  oscillator  strength  of  the  C"-A'  transition  relative  to  the  B-A 
transition.  As  mentioned  in  the  introduction,  the  calculated  radiative 
lifetime  for  the  C"-A'  transition,  4.3  us,  conflicts  with  some  of  the  data  on 
quenching  of  C"-state  emission.  The  quenching  data  suggests  that  the  lifetime 
of  the  C"-state  is  closer  to  40  ns.  This  lifetime  would  imply  a  much  stronger 
absorption  oscillator  strength  for  the  transition.  The  discrepancy  will  be 
addressed  further  in  the  discussion  session;  however,  the  primary  effect  of  an 
increased  oscillator  strength  for  the  C"-A'  transition  would  be  a  lowering  of 
the  detection  limit  for  the  A' -state  by  approximately  a  factor  of  one  hundred. 
Thus,  our  measurements  would  give  a  result  which  is  an  upper  limit  on  the 
A' -state  concentration  and  hence  a  lower  bound  on  the  effective  loss  rate. 

An  LIF  spectrum  of  the  region  from  700  to  710  nm  is  shown  in  Fig.  37.  The  2-0 
and  3-1  bands  of  the  C"-A'  transition  should  be  excited  in  this  wavelength 
region.  Though  absorption  features  are  observed  in  this  region,  they  can  be 
assigned  to  the  N2<B-A,  8-6)  transition.  This  assignment  is  confirmed  by 
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Figure  37.  LIF  excitation  spectrum  for  the  region  from  700  to  710  nm.  Features 
of  the  B-A,  8,6  and  7,5  transition  are  identified.  Arrows  show  the 
expected  location  of  the  HIR  3,1  and  2,0  bands.  The  excitation 
point  for  the  spectrum  of  Fig.  38  is  shown. 

resolving  the  fluorescence  excited  by  the  laser.  A  typical  resolved  emission 
spectrum  is  shown  in  Fig.  38.  The  only  discernible  emission  features  are  due 
to  transitions  originating  from  v=8  and  v=7  of  the  B-state.  Population  of  N£ 
(B,  v=7)  is  due  to  vibrational  relaxation  of  B,  v=8  by  the  Ar  buffer.  Thus,  no 
emission  from  the  ^(C'-A')  system  is  excited  by  the  laser  within  the  detection 
limits  of  the  system.  The  detection  sensitivity  used  in  the  LIF  spectrum  of 
the  700  to  710  nm  region  is  a  factor  of  1000  better  than  than  the  sensitivity 
used  in  recording  the  spectrum  of  the  N2(B-A,  3-0)  band.  Observation  of  the 
poorly  populated  and  weakly  absorbing  v=6  level  of  the  A-state  illustrates  the 
detection  limits  of  the  system.  Comparison  of  the  noise  levels  observed  in 
this  spectrum  with  the  signal  levels  obta  d  in  the  spectrum  of  Fig.  37, 
corrected  for  the  differences  in  the  absorption  strengths  for  the  B-A,  3-0  and 
C"-A' ,  2-0  bands,  establishes  a  detection  sensitivity  of  approximately 
1  a  10**  molecules  cm--*. 
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Wavelength  (nm) 

Figure  38.  Resolved  LIF  spectrum  for  one  of  the  absorption  bands  of  Fig.  37. 

The  features  present  can  all  be  identified  as  originating  from  the 
B-A  transition.  Excitation  of  the  8,6  band  is  clearly  indicated. 
Emission  from  B,  v=7  is  probably  due  to  vibration  relaxation  of  the 
B-state. 

The  simple  numerical  model  used  to  predict  the  expected  N2(A' )  densities  at  the 
observation  point  in  the  discharge  flow  reactor  can  be  used  to  determine  a 
lover  limit  on  loss  rate  (kg)  for  N2<A')  consistent  with  a  concentration  at 
our  detection  limit.  The  traces  in  Fig.  39  show  the  expected  temporal  profile 
for  N2(A')  when  only  wall  loss  (Trace  a)  and  a  loss  rate  of  3200  s"^  (Trace  b) 
is  included  in  the  model.  When  the  loss  rate  of  3200  s-^  is  incorporated  in 
the  model  an  N2(A')  concentration  equivalent  to  the  estimated  detection  limit 
is  predicted  at  the  detection  point.  Hence,  we  can  place  an  upper  limit  on  the 
A' -state  lifetime  of  approximately  0.3  ms.  In  our  experiments  we  explored  a 
range  of  N2/Ar  mixture  ratios  (1  to  20  percent  N2),  discharge  conditions  (300 
to  600  VDC) ,  and  pressures  (2  to  5  Torr).  In  addition,  we  passed  N2/Ar  and 
N2/Ar/He  mixtures  directly  through  the  DC  discharge.  This  method  of  production 
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Figure  39.  Model  predictions  for  the  temporal  evolution  of  the  ^(A' )  popula¬ 
tion  in  the  discharge  flow  reactor.  Trace  a  shows  the  behavior 
expected  if  wall  loss  was  the  only  mechanism  responsible  for 
quenching  of  the  A' -state.  Trace  b  shows  the  behavior  expected  for 
a  loss  rate  of  3200  s“l.  This  rate  produces  a  concentration  of  A' 
at  our  detection  limits. 

enhances  the  observed  ^(C-A')  emission  intensity  at  the  expense  of  producing 
N-atoms  and  perhaps  other  metastable  states  of  N2-  In  addition,  He/Ar  mixtures 
were  co-discharged  and  N2  added  downstream.  Under  no  set  of  conditions  was  LIF 
on  the  C"-A'  transition  observed. 

2.2.5  Discussion 

Ve  have  conducted  a  series  of  experiments  to  determine  the  gas-phase  lifetime 
of  the  A' -state  of  N2.  Though  the  LIF  calibration  method  employed  in  the  study 
was  not  precise,  the  experiments  do  infer  a  lifetime  for  the  -‘L  state  which  is 
much  shorter  than  expected  for  a  spin-forbidden  decay  process.  Our  experiments 
do  not  provide  any  direct  evidence  for  any  of  the  possible  decay  mechanisms. 
However,  some  indirect  evidence  provides  support  for  a  collisional  decay 
channel  through  the  B^ITg  state. 
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An  expanded  potential  energy  diagram  for  the  A' -state  and  near  isoenergetic 
levels  of  the  B^II  and  A^I  state  is  shown  in  Fig.  40.  The  diagram  shows  that 
the  lowest  vibrational  level  of  the  A' -state  is  near-resonant  with  v=10  of  the 
Ilg  state.  Also,  the  barrier  to  predissociation  for  the  A'-state  occurs  at 
v'=13  of  the  B-state.  The  latter  observation  provides  evidence  for  invoking 
curve  crossing  to  the  A'-state  as  the  mechanism  for  predissociation  of  the 
B-state.  Thus  both  collisional  and  spin-orbit  coupling  to  the  B-state  might  be 
evident  in  the  spectrum  of  levels  10-12.  A  recent  study  of  the  high-resolution 
LIF  spectrum  from  the  ^(B^  Ilg-A^  Zu+ ,  10-6)  transition  revealed  no  apparent 
perturbations  in  the  rotational  levels  of  the  B-state. 51  In  addition,  no  other 
studies  have  reported  any  perturbations  in  the  spectra  of  levels  10-12.52 
Thus,  there  appears  to  be  no  direct  evidence  from  the  spectroscopy  of  the 
B-state  to  indicate  spin-orbit  coupling  to  the  A'-state. 


Figure  40.  Expanded  energy  level  diagram  of  the  A' state  redrawn  from  the 
results  of  Partridge  et  al.56  The  figure  illustrates  the 
potential  for  coupling  of  the  A'  and  B  states. 
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The  evidence  for  collisional  coupling  of  the  B-state  to  the  A'-state  come  from 
measurements  of  the  energy  pooling  rate  coefficient  to  form  the  B-state.  If 
the  energy  pooling  reaction  is  studied  using  an  Ar  buffer  gas  much  of  the 
information  on  the  distribution  of  vibrational  levels  populated  in  the  B-state 
is  washed-out  by  rapid  vibrational  relaxation.  If  He  is  employed  as  a  buffer 
gas  vibrational  relaxation  is  less  efficient  and  much  of  the  nascent  vibra¬ 
tional  distribution  is  preserved  in  the  emission  spectrum.  A  typical  spectrum 
of  emission  from  the  energy  pooling  reaction  from  560  to  720  nm  using  He  as  a 
buffer  gas  is  shown  in  Fig.  41.  In  contrast  to  the  spectrum  of  Fig.  38,  a 
strong  enhancement  in  the  population  of  v=10  of  the  B-state  is  observed.  This 
enhancement  has  been  previously  observed. ^3  Thus,  it  would  appear  that 
collisional  coupling  from  the  A'-state  to  the  B-state  may  be  a  mechanism  for 
the  decay  of  the  A'-state. 


WAVELENGTH  (nm) 

B-0758* 


Figure  41.  Visible  emission  spectrum  of  N2  due  to  N£(A)  energy  pooling  in  a 
dominantly  He  buffer  gas.  The  enhancement  of  the  population  of 
v=10  of  the  B-state  is  indicated  in  the  inset  plot  of  the 
vibrational  distribution. 
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Additional  evidence  for  the  coupling  comes  from  the  magnitude  of  the  rate  coef¬ 
ficients  to  form  the  B-state  and  the  C"-state  from  the  energy  pooling  reaction. 
The  energy  pooling  rate  coefficient  for  the  direct  production  of  N2(B)  from  two 
N2(A,  v=0)  molecules,  as  measured  by  Piper, ^3  js  7.7  +  1.1  x  10“H  cm3 
molecule-!  s-!.  This  rate  coefficient  can  be  compared  with  his  measurement  for 
the  formation  of  the  C"-state,  8.1  +  2.3  x  10-!!  cm3  molecule-!  s-!.  Thus,  it 
is  possible  that  a  significant  fraction  of  the  direct  production  of  B-state  is 
actually  derived  from  production  of  the  C"-state  followed  by  radiative  cascade 
to  the  A' -state  and  collisional  coupling  to  the  higher  B-state  vibrational 
levels.  However,  some  direct  production  of  the  B-state  must  occur.  Piper's 
observation  that  energy  pooling  of  vibrationally  excited  N2(A)  favors  pro¬ 
duction  of  lower  B-state  levels  is  not  consistent  with  the  cascade  mechanism. 

In  addition,  pooling  of  N2(A,  v=0)  molecules  also  shows  enhanced  production  of 
B,  v= 2. 

The  upper  limit  on  the  stability  of  the  A'-state  is  uncertain,  to  a  large 
degree,  due  to  the  apparent  discrepancy  in  the  lifetime  of  the  C"-state. 
Predissociation  of  the  C"-state  can  be  excluded  as  the  cause  for  the  apparent 
reduced  lifetime.  The  energy  pooling  rate  coefficient  would  have  to  be 
significantly  larger  than  gas  kinetic  in  order  to  produce  a  C"-state 
population  for  which  the  predissociation  rate  was  100  times  greater  than  the 
fluorescence  rate.  Similarly,  direct  coupling  of  the  C  and  C"  states  is 
unlikely.  At  the  pressures  employed  in  most  of  the  experiments  which  have 
investigated  this  system  radiative  decay  of  the  C-state  should  be  much  more 
efficient  than  collisional  energy  transfer.  Further  investigation  of  this 
problem  is  warranted. 

Our  experiments  in  which  we  vary  the  N2  concentration  indicate  that  N2  is  not 
primarily  responsible  for  coupling  of  the  two  manifolds.  At  the  lowest  N2 
concentrations  employed  a  gas-kinetic  coupling  rate  coefficient  would  be 
required  to  couple  the  levels  and  significant  variations  in  the  B-state,  v=10 
enhancement  would  have  been  observed.  It  would  appear  more  likely  that  Ar  or 
He,  which  are  both  present  in  much  greater  concentrations,  may  be  responsible 
for  the  coupling.  The  production  of  ^(A')  in  either  a  supersonic  jet  or  in  a 
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rare  gas  matrix  may  reduce  the  collisional  coupling  to  the  B-state  and  enable  a 
more  rigorous  determination  of  the  lifetime  of  this  state. 

2.2.6  Conclusions 


We  have  observed  production  of  the  A'^Zg+  state  of  N2  from  the  energy  pooling  of 
N2(A^Eu+)  molecules.  Laser-induced  fluorescence  on  the  C^Ily-A' 5 transition 
failed  to  detect  the  presence  of  A' -state  molecules.  Based  on  simple  kinetic 
arguments,  we  have  been  able  to  place  an  upper  limit  on  the  lifetime  of  the 
A'-state  under  our  experimental  conditions  of  0.3  ms.  Observation  of  enhanced 
emission  from  v=10  of  the  B^llg  state  suggests  that  collisional  coupling  of  the 
A'  to  the  B-state  is  responsible  for  the  decay  of  this  state. 

The  production  and  storage  of  the  A'-state  under  cryogenic  conditions  via  matrix 
isolation  may  enhance  its  stability.  The  energy  gap  between  v=10  of  the  B^IIg 
state  and  v=0  of  the  k'^Zg+  state  is  approximately  650  cm~l  which  is  much 
smaller  than  kT  at  4K.  Thus  collisional  coupling  at  this  temperature  should  be 
quite  inefficient.  In  addition,  the  magnitude  of  spin-orbit  coupling  to 
adjacent  electronic  states  is  also  reduced  at  low  temperatures.  The  selection 
rules  for  spin-orbit  coupling  to  adjacent  electronic  states  is  also  reduced  at 
low  temperatures.  The  selection  rules  for  spin-orbit  coupling  require  that 
total  angular  momentum  be  conserved  (AJ=0)  and  that  spin  and  orbital  angular 
momentum  have  the  selection  rules  (Reference  54): 

AL  =  0,  ±  1  AS  =  0,  +  1  (39) 

The  total  angular  momentum,  J,  is  the  sum  of  the  spin,  orbital,  and  rotational 
contributions  (neglecting  nuclear  spin).  At  cryogenic  temperatures  (and  in  the 
matrix)  rotational  motion  is  frozen  into  the  lowest  rotational  state  and 
rotational  transitions  are  not  possible.  Hence,  the  spin-orbit  coupling 
selection  rules  require  that  AS  =  AL.  Since  there  are  no  other  nearby  quintet 
states,  spin-orbit  coupling  of  the  5E-state  requires  AS  =  +1.  Hence,  the  ^Z 
state  can  only  couple  to  or  states.  Thus,  the  B^IIg  state  is  the  only 
state  which  is  allowed  couple  to  the  A'5£g+  state.  The  reduced  degree  of 
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rotational  motion  imposed  by  the  matrix  results  in  a  reduced  number  of  pathways 
for  spin-orbit  interaction  and  hence  promotes  the  stability  of  the  state. 

Thus,  the  challenge  to  the  further  stabilization  of  the  A'^E  state  is  to  form 
and  stabilize  this  molecule  in  a  cryogenic  matrix.  In  the  paragraphs  below  we 
discuss  the  previous  work,  on  cryogenic  trapping  of  excited  N2  and  the  evidence 
which  suggests  that  the  A' -state  can  be  matrix  stabilized  in  a  form  useful  for 
propulsion. 

The  concept  of  energy  storage  by  cryogenic  trapping  of  energetic  radical  and 
excited  state  species  in  matricies  arose  some  30  years  ago  and  was  vigorously 
pursued  in  an  intense  research  program  at  the  National  Bureau  of  Standards  (now 
NIST). ^5, 56  xhe  program  focussed  on  the  trapping  of  ground  state  atomic  and 
molecular  radical  species.  However,  in  the  course  of  the  research,  several 
investigators  observed  complex  fluorescence/phosphorescence  effects  which  were 
ascribed  to  cryogenic  trapping  and  subsequent  energy  transfer  involving  excited 
state  species.  Much  of  this  early  work  can  be  found  in  the  papers  of  Broida  and 
coworkers  (Refs.  57-60  and  citations  therein)  and  Brocklehurst  &  Pimentel. 61  in 
brief,  these  investigators  introduced  flowing  discharged  N2/rare  gas  mixtures 
into  a  1  quid  He  cooled  trap  where  the  discharged  species  were  condensed.  Upon 
cessation  of  the  discharge,  the  condensed  samples  were  observed  to  glow  in  the 
visible,  with  characteristic  decay  times  on  the  order  of  seconds  to  minutes, 
indicate  lg  rather  lengthy  storage  of  the  trapped  electronic  energy.  The  solid 
afterglow  was  also  observed  to  follow  a  rather  complex  temperature  dependence 
upon  warning  the  sample  from  A. 2  K  to  35  K,  indicating  that  the  energy  decay 
pathways  were  a  function  of  matrix  material  and  matrix  diffusion  times. 

More  recently,  the  excitation  of  emissions  in  N2  at  cryogenic  temperatures  was 
summarized  and  further  addressed  by  Dressier  and  coworkers. 62-64  They  showed 
that  these  emissions  were,  for  the  most  part,  independent  of  the  excitation 
source  employed.  Furthermore,  they  were  all  related  to  the  presence  of  the 
excited  and  states  of  nitrogen  atoms.  The  band  systems  observed  can  be 
assigned  to  the  following  processes:  a)  phonon  assisted  emission  on  the  ^D-^S 
and  ^P-^D  transitions  with  characteristic  lifetimes  of  37s  and  approximatly  1ms 


67 


respectively,  b)  a  S tokes-shi f ted  transition  corresponding  to  emission  on  the 
2D-^S  transition  coupled  to  excitation  of  a  2300  cm-^  vibron  mode  in  the  N2 
crystal  lattice,  and  c)  Stokes  and  anti-Stokes  sequence  band  emission 
corresponding  to  the  formation  and  radiation  of  N(2P)  and  N(2D)  coupled  with 
high  vibrational  transitions  in  the  N2  lattice.  These  features  are  summarized 
in  Fig.  42  for  both  the  2D-4S  and  2P-2D  transitions. 

The  appearence  of  these  vibrational  satellites  has  been  ascribed  to  the  quench¬ 
ing  of  matrix  created  N2(A2£U+)  by  N-atoms  to  produce  either  N(2D)  or  N(2P), 
with  the  remaining  energy  taken  up  in  the  vibrational  modes  of  the  matrix. 

Vibron  levels  as  high  v=14  may  be  excited  in  the  formation  of  N(2D)  and  v=8  in 
the  formation  of  N(2P).  These  limits  are  consistent  with  the  cut-off  in 
vibrational  levels  observed  corresponding  to  excitation  of  the  respective 
states.  Emission  from  the  N2<A-X)  Vegard-Kaplan  bands  was  readily  observed  by 
Vegard  in  excited  N2~doped  solid  rare  gases, 65  but  was  observed  to  weaken  if  the 
matrix  was  irradiated  with  electrons.  It  was  hypothesized  that  the  gradual 


Figure  42.  Energy  level  diagram  and  transitions,  of  the  visible  and 
near-infrared  luminescence  of  N-atom-doped  solid  N2. 
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increase  in  N-atom  density  accompanying  the  irradiation  resulted  in  the 
well-known  quenching  of  N2(A)  by  N2.66-68  Thus,  the  lifetime  of  the  excited 
triplet  appears  to  be  directly  related  to  the  N-atom  density  in  the  matrix. 

The  lifetimes  reported  for  N(^D)  and  N(3p)  in  solid  N2  represent  a  reduction  of 
~103  and  ~104  over  the  respective  gas-phase  lifetimes. ^9  Dressier  and 
coworkers  attribute  the  lifetime  shortening  to  a  matrix-induced  increase  in  the 
transition  rate.  However,  no  evidence  would  appear  to  support  this  interpreta¬ 
tion  over  a  collisional  matrix  relaxation  process.  This  distinction  is  critical 
since  radiation  trapping  can  be  used  to  retain  stored  energy  in  a  thick  matrix 
if  an  increase  in  transition  rate  is  responsible  for  the  shortened  lifetimes. 
Otherwise,  this  excitation  is  lost  to  the  matrix  heat  bath. 

The  lifetime  of  the  N(3D-4S)  transition  and  the  N2(A-X)  bands  as  a  function  of 
matrix  material  and  temperature  were  studied  by  Tinti  and  Robinson. The 
lifetime  of  the  3D-4S  transition  appeared  to  lengthen  as  the  rare  gas  matrix 
partner  became  lighter.  Lifetimes  of  340  s,  18  s,  and  7  s  were  measured  for  Ne, 
Ar,  and  Kr  matricies  respectively.  The  lifetime  of  N2(A,  v=0)  showed  similar 
lengthening  with  rare  gas  partner.  However,  unlike  the  N-atom  transitions,  the 
N2<A-X)  transition  lifetime  was  nearly  identical  to  the  reported  gas  phase 
lifetime,  within  the  uncertainties  of  the  respective  measurements.  This  can  be 
attributed  to  the  reduced  probability  of  spin-orbit  coupling  between  the 
molecule  and  the  host  matrix.  This  is  an  extremely  important  observation.  The 
magnitude  of  the  matrix  elements  which  control  these  couplings  are  directly 
related  the  the  viability  of  the  N2(A'3Eg+)  state  as  an  energy  storage  medium. 
The  data  of  Tinti  and  Robinson  also  show  that  the  lifetime  is  independent  of  N2 
dosing  levels  up  to  1  percent.  No  data  was  reported  for  higher  dosing  levels. 
The  data  also  agrees  well  with  the  calculations  of  Happer  and  coworkers^3  who 
showed  that  the  lifetime  of  pure  spin-polarized  N2  crystals  is  quite  short  due 
to  the  magnetic  dipole-dipole  interactions  in  the  crystal.  However,  they  also 
noted  that  the  stability  of  high-spin  states  may  be  enhanced  by  dilution  in  a 
diamagnetic  matrix  like  solid  hydrogen.  Both  theory  and  experiment  point  to  the 
dilution  in  a  light  atom  matrix  as  a  means  of  stabilizing  these  high  spin 
systems.  Therefore,  we  believe  that  this  system  merits  further  investigation. 
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